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INFORMATION 


D-N-L#17 


Two letters from Halvor DEVOLD, Norway: 

(with some good reasons to be printed on the information page!!) 

Dear Josef, ... Thank you for the exe-file of my KEPLER-movie-picture and please forgive me for this 
delayed response! Time has simply run away for me since I came back from Plymouth last summer. 
Our new school reform, Reform94, introduced the graphic calculator as a compulsory tool, to be used 
on all first-year courses in mathematics. Most of our math teachers were taken by surprise, and re- 
acted accordingly, that is in the negative ---. 

Hence, I have had - and still have - a hard time in persuading colleagues to see the advantage of such 
tools, - and to teach my nearest colleagues how to use graphic colleagues as well as DERIVE. We 
are just a few math-teachers in this country - it seems - who have realized the new situation - and the 
new opportunities - that has arisen in maths teaching, and who are also in regular contact - through a 
BBS. (Firda BBS) 

In order to improve upon this situation I have borrowed your idea and made a printed bulletin with 
extracts from the same BBS-conference. (Its name is DRIVE) 

You will find the first issue of this printed bulletin attached. It is of course by no means meant to com- 
pete with your bulletin of the DUG, but solely directed to a Norwegian public, and it is not - as you will 
see - devoted to DERIVE alone, but meant to be informative on both graphic calculators and DERIVE. 
In fact, with DERIVE 3, the gap between graphic calculators and CAS seems to become narrower, - 
that is in their user interface. Hence it will be easier to persuade people to use both, and to see and 
judge for themselves when to use DERIVE, and when a graphic calculator is sufficient. Both have their 
virtues Best regards. 

More information about FIRDA & DRIVE: 

Halvor Devoid, Maneweien 10, N-4900 Tvedestrand, Norway 

Dear Josef, I am now working on a project that includes the use of AcroSpin: 

I bought AcroSpin from MathWare, US recently in order to explore its potential as a means to extend 
the visual power of DERIVE. My first aim was to see if it could be used to present the five regular poly- 
hedrons in stereo vision 


DNL. I think, that some of you have wondered about the Option to save a 3D-plot in AcroSpin format. Here you 
can find the answer. I can assure you, Halvor produced a wonderful file. I asked him to inform us how to obtain 
AcroSpin. Halvor answered immediately: 

You can order AcroSpin at MathWare, P.O.Box 3025, Urbana, IL, 61801, (800) 255 2468, 

FAX (217)384 7043 

The price is US$ 30.00. MathWare offers a second useful program called GyroGraphics (US$ 89.00). 
GyroGraphics supports AcroSpin. 

Denken Sie bitte an das 


1 • Deutsche DERIVE User Group Treffen 
im Rahmen der MNU 1995 
Niirnberg, 9. - 13. April 1 995 

Das Treffen findet am Mittwoch, den 12.April 1995 
um 14 Uhr staff. 

Die Raumlichkeit erfahren Sie im Tagungsbiiro. 

Fakultat der Wirtschafts- und Sozialwissenschaften, 
Uni Erlangen-Nurnberg, Lange Gasse 20 


Ich wiirde mich freuen, viele von Ihnen wieder zu sehen, bzw. kennen zu lernen. 
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LETTER OF THE EDITOR 


Liebe DUG-Mitglieder, 

Ich begriifie Sie alle recht herzlich zu Beginn des 
5. DUG-Jahres. Ganz besonders herzlich willkom- 
men heifie ich die neuen DUG-Mitglieder. So ha- 
ben wir seit kurzem auch eine „diplomatische” 
Vertretung in Neuseeland und in Kolumbien. 

Ich habe einige Anfragen bezuglich KEPLER.EXE 
erhalten. Leider ist GRASP - mit diesem Prasenta- 
tions- und Animationsprogramm wurde KEPLER 
erzeugt - nicht mit alien Grafikkarten problemlos 
kompatibel. Oft hilft es, ein wenig mit den Video- 
modi zu experimentieren. H.Devold, der Schopfer 
von KEPLER, hat mit einem Einsatz von AcroSpin 
- DERIVE3 Anwender kennen AcroSpin zumindest 
von Transfer-Save-AcroSpin im 3D-Menu - wieder 
eine wahrhaft bezaubernde Show geschaffen. 

Diese Ausgabe des DNL ist zu einem betrachtlichen 
Teil dem DERIVE - Einsatz in osterreichischen 
Schulen gewidmet. Ich bedanke mich bei alien 
Kollegen fur ihre Bereitschaft, uns einen Blick in 
ihren Schulalltag werfen zu lassen und auch fur 
ihre Geduld bei meinen lastigen Nachfragen. Allen 
"glucklichen?” Nichtlehrern unter Ihnen verspre- 
che ich, dass die nachsten Ausgaben nicht nach 
Kreide und Tafel riechen werden. Ich habe fur 
dieses Exemplar einen Index fur alle bisherigen 
Ausgaben des DNL vorbereitet, aber wie so oft, der 
Platz reichte trotz des deutlich gestiegenen Um- 
fangs wieder nicht aus. 

Ich hoffe, in nachster Zukunft auch uber e-mail 
erreichbar zu sein. Die Ausziige aus dem BBS und 
aus der derive-news@mailbase.ac.uk. in Birming- 
ham haben mir Appetit gemacht. Eine digitale 
Derive-Konversation - Lehrstunde konnen Sie ab 
Seite 7 verfolgen. 

Im nachsten DNL will ich anstelle des Beitrags von 
Thomas Weth, dem ich an dieser Stelle ganz beson- 
ders fur seine zuverldssigen Zusendungen danken 
mochte, einen Artikel von Peter Baum, Kassel, 
veroffentlichen, der aus seiner Beschaftigung mit 
Thomas Weths „Kurvenlexikon" entstanden ist. Es 
freut einen Herausgeber wirklich, wenn sein Bulle- 
tin nicht nur abonniert, sondern auch tatsachlich 
gelesen wird und wie ersichtlich, sogar Anregun- 
gen fur eigene Arbeiten liefert. 

Abschliefiend mochte ich Sie noch auf die MNU in 
Niirnberg - mit dem deutschen Usergrouptreffen - 
und auf die DERIVE DAYS Diisseldorf hinweisen. 

Mit vielen Griifien bis zum Sommer ~ ^ — 5 

* ^ 

(PS.: Ihre Mitgliedsnummer finden Sie auf dem 

A dressaufkleber) 


Dear DUG-Members, 

My best regards to you all at the beginning of 
the 5th DUG year. My special welcome to the 
new DUG members. So we have opened a 
"diplomatic" representation in New Zealand 
and in Colombia. 

I received some questions concerning 
KEPLER.EXE. Unfortunately it seems that 
GRASP - the presentation- and animation 
program used to produce KEPLER - is not 
compatible with some graphics adapters. It 
might help to experiment with the video modes 
of your graphics card. H.Devold, the creator of 
KEPLER, supported by AcroSpin - DERIVE 3 
users know AcroSpin at least from Transfer- 
Save-AcroSpin in the 3D Menu - produced a 
new charming show. 

A big part of this issue is dedicated to the Aus- 
trian DERIVE scene. I want to thank all my 
colleagues for their willingness to give us in- 
teresting insights in their daily life at school. 
And I thank them for their patience when I had 
many boring questions about their contribu- 
tions. But I promise all the "fortunate?" non- 
teachers among you, that the next DNL will not 
smell of chalk and blackboard. For this issue I 
have prepared an index for all the DNLs from 
DNL#1 - DNL#16, but I ran out of paper, de- 
spite the now 42 pages contents. 

I hope to be available by e-mail in the near 
future. The excerpts from the BBS and from 
the derive-news@mailbase.ac.uk have given 
me an appetite. You can follow a digital 
DERIVE conversation - lecture on page 5. 

Instead of the next chapter of Thomas Weth 's 
"Lexicon of Curves" - let me thank him for his 
reliability in delivering his contributions - 1 want 
to publish an article by Peter Baum, Kassel, 
which emerged from dealings with Thomas 
Weth's "Lexicon". Its really a pleasure for an 
editor, to see that people do not only renew 
their subscriptions but actually read the journal 
and find in it ideas for their own investigations. 

With my best wishes until summer 

Sincerely 



(You can find your DUG membership number 
on the address label.) 
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The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE User Group. It is pub- 
lished at least four times a year with a con- 
tents of 40 pages minimum. The goals of 
the D-N-L are to enable the exchange of 
experiences made with DERIVE as well as 
to create a group to discuss the possi- 
bilities of new methodical and didactical 
manners in teaching mathematics. 


Preview: Contributions for the next issues 

Stability od systems of ODEs, Kozubik, SLO 
Prime Iterating Number Generators, Wild, UK 

Graphic Integration, Probability Theory, Linear Programming, Bohm, AUS 

Continued Fractions and the Bessel Functions, Cordoba a.o., ESP 

LOGO in DERIVE, Lechner, AUS 

DREIECK.MTH, Wadsack, AUS 

IMP-Logo and Misguided Missiles, Sawada, HAWAII 

"Reverse Discussion" of Curves, Reichel, AUS 

Reichel - Klingen - Bohm - Splines, a triathlon, AUS & GER 

3D-Geometry, Reichel, AUS 

Parallel- and Central Projection, Bohm, AUS 

Conic Sections, Fuchs, AUS 

Muller's Method to solve univariate equations, Speck, NZL 
Vector and Vector Indices Sorting, Biryukow, RUS 
Algebra at A-Level, Goldstein, UK 
Tilgung fremderregter Schwingungen, Klingen, GER 
A Utility file for complex dynamic systems, Lechner, AUS 
and 

Setif, FRA; Vermeylen, Belgium; Lymer, FRA; Leinbach, USA; Baum, GER; 
Weth, GER; Wiesenbauer, AUS; Keunecke, GER; 

and messages from the derive-news@mailbase.ac.uk 


Editor: Mag. Josef Bo hm 
A-3042 Wiirmla 
D'Lust 1 
Austria 

Phone/FAX: 43-(0)2275/8207 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the D-N-L. 
It must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
D-N-L. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE Newsletter will be. 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wiirmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
Herstellung: Selbstverlag 
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R.Larham, Abbots Langley, UK 

Alfonso J.Poblacion is mistaken about the curve (DNL#16, p6) given by the parametric equations: 

[sin(0.99t) - 0.7 cos(3.01t) , cos(l.Olt) + 0.1 sin(15.03t)] 

It is clearly closed as it is continuous and both the x and y-coordinates are periodic with a common 
period of 200tt. To demonstrate this try plotting the curve for t e [20071, 200 tt+27t] over the plot for 
t e [0,27i]. Alternatively try simplifying: 

[sin(0.99t) - 0 . 7cos (3 . Olt) , cos (1 . Olt) + 0 . lsin (15 . 03t) ] - - 
[sin(0.99(t + 20071)) - 0.7cos(t + 200 ti) ) , cos (1 . 01 (t + 20071)) + 

+ 0 . lsin (15.03 (t + 200 tt) ) ] 

If this curve does not appear to close when plotted for a large enough range for t , it is an artefact of the 
plotting algorithms and not of the curve itself. 

In fact it is easy to see that any curve which can be represented in parametric form by: 

[A sin (at) - B cos (Pt) , C cos (yt) + D sin (St)] 
is closed, when a, (3, y, 5 are rational multiples of one another. 

Dr. A.Lvubimov, St- Petersburg, Russia 

Dear Sir, thank you very much for your kind letter, for invitation to join DUG and especially for DNL 
you have sent me. It is a very useful and fascinating issue. 

I have met some difficulties with the sending of my membership due (this is the reason of my delay 
with the answer) because our young banks have the weak connections with Europeans ones yet. So I 
had to transfer my due on 1995 with my friend who was going abroad. He had to send my due and the 
application form to you last week from London. Please let me know if you did not receive them. 

I gladly met your information about Sergey Biryukow from Moscow. Please send me his mail ad- 
dress and/or e-mail if possible so that I could contact him. 

I am very interested in cooperation with you, too and I also know the strong background of Euro- 
pean mathematics and teachers. Some ideas of my High School course I wrote are taken from the bril- 
liant "Elementarmathematik vom Hoheren Standpunkte aus, l.Band" by Felix Klein that had been 
translated into Russian just ten years ago. In addition I am going to study German. 

I hope our cooperation will be useful and fruitful. Best wishes for Christmas and New Year. 
A.Lyubimov 

A.Kropf, Hildesheim, Germany 

Message 3989: From AL POOR to AL RICH about ITERATES: 

Authorf(x) := x A 2 and ITERATES (f(x) ,x,x, 3). Simplifying yields [ x, x 2 , x 8 , x 127 ] . Is this 
a known bug - or what's the matter? Sincerely 
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DNL: I don't know, why you have x 127 


, but your result made me curious. Compare the two listings: 


DERIVE Version 3.01 

u : =x A 2 

ITERATES (u, x, x, 3) 

[x, x A 2 , x A 4 , x A 8 ] 

F (x) : =x A 2 

ITERATES (F (x) , x, x, 3) 

[x, x A 2, x A 4, x A 8] 
v : =2*x+3 

ITERATES (v, x, x, 3) 

[x, 2*x+3, 4*x+9, 8*x+21 ] 
G (x) : =2 *x+3 
ITERATES (G(x) , x, x, 3 ) 

[x, 2*x+3, 4*x+9, 8*x+21 ] 
ITERATES (v, x, 2, 3) 
[2,7,17,37] 

ITERATES (G(x) , x, 2 , 3 ) 
[2,7,17,37] 


DERIVE Version 3.04 

u : =x A 2 

ITERATES (u, x, x, 3) 

[x, x A 2, x A 4, x A 8] 

F (x) : =x A 2 

ITERATES (F (x) , x, x, 3) 
[x,x A 2,x A 8,x A 128] 
v: =2*x+3 

ITERATES (v, x, x, 3) 

[x, 2*x+3, 4*x+9, 8*x+21 ] 

G (x) : =2*x+3 
ITERATES (G (x) , x, x, 3) 

[x, 2*x+3, 8*x+21, 128*x+381] 
ITERATES (v,x, 2, 3) 
[2,7,17,37] 

ITERATES (G (x) , x, 2, 3) 
[2,7,17,37] 


I don't know the reason, why the "software smithies" have changed the ITERATES-function, but ob- 
viously there is now a difference whether you enter an expression - say v = 2x+3 or a function 
G(x) = 2x+3 as iteration rule. 

In the first case we express x„+i = 2 x„ + 3 with xO = x, in the second case we do as follows: 


xO = x 

xl = G(x0) = 2x + 3 

x2 = 2 G(xl) + 3 = 2 (2 xl + 3) + 3 = 2 (2 (2x + 3) + 3) + 3 = 8x + 21 

x3 = 8 G(x2) + 21 = 8 G(2 x2 + 3) + 21 = 8 (2 (8x + 21) + 3) + 21 = 128x + 381. 


I'll try to find an answer from A1 Rich for you (and for us!). Josef. 

(Mr Kropf, please give me your address, that I can send you an answer.) 

DERIVE 6 makes no difference and gives the results as version 3.01 did a couple of years 
ago. 

There was another problem occurring with the ITERATES-function. You know Johannes 
Wiesenbauer, our sophisticated specialist in number theory. In December he called me very disap- 
pointed because his beloved POWERMOD-function from DNL#14 didn't run in DERIVE 3.00, al- 
though he had no problems with the Beta-version. I contacted SWHH and A1 Rich faxed an answer: 


Dear Josef, thank for reporting the problem DERIVE and DERIVE XM Version 3.00 have 
simplifying calls on the INV function. I am happy to report that the problem has been found 
and is now corrected by Version 3.01. 

I will request that SWHE send a Version 3.01 diskette to Dr. Wiesenbauer so he can verify 
that all his functions now work correctly. 

I sincerely apologize for any inconvenience this bug may have caused you and the other 
members of the DERIVE User Group. 


D 
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As we change DERIVE to make it more powerful and efficient, new bugs are inevitably intro- 
duced despite our best efforts. Before releasing a new version, I run it through a test suite of 
over 7,000 math problems. Since I add to the test suit problems designed to detect the bugs 
you and others report (e.g. the INV function will now be added), at least the same bug should 
not occur again! (You might put a note to this effect in the DERIVE Newsletter.) 

What I've done now, and many thanks, At for giving us an insight into your work. 


M.Th.Lanqsch, Wittlich, Germany 


Gleichungen 3. Grades scheinen in der Luft zu liegen. Im Schulbuch fur die II.KIasse steht 


die Aufgabe, ein s zu finden, so dass fur alle x > s gilt: 


1 


-x— 1 


< 


1 


1000 


, das fuhrt auf 


x 3 - x - 1001 = 0. Das Losungsbuch gab ohne Kommentar s = 10 an (!).lch anderte den Nen- 
ner in x 3 - x + 10. 


Nun waren wir bei x 3 - x - 990 = (x - 10) (x 2 + lOx - 99) = 0 und kannten die einzige reelle 
Losung. Doch die Schuler wollten wissen, wie man Gleichungen 3. Grades lost. Mit den Car- 
danischen Formeln kamen wirzu: 


x = 3 . 495 + , 495 2 


v3y 


+ 3 ,495- J495 2 ■ 


1. 

v3y 


, doch nun kamen wir mit Termumformun- 


gen nicht weiter. DERIVE macht daraus mit approX zwar 10, mit Simplify aber nur: 


x= 3 495+ 


299V222 


+ 3 495- 


299V222 


- warum? 


Zur Erganzung sei erwahnt, dass 495 + 


299V222 


5 


ist, und damit trivialerweise 



Mit freundlichen GruBen 

Recent versions of DERIVE have no problems simplifying the nested roots! 
3 

#1: SOLVECx - x - 990, x) 

#2: x = -5 - J74-L v x = -5 + J74-L v x = 10 


495 + 4 

2 

495 - 

' 1 ' 

3'' 

'1/3 

+ 

f 

495 - 4 

2 

495 - 

' 1 ' 

3Y 

\ 


, 3 , 

J . 




, 3 > 

// 


#4: 


10 
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The Tl-Calculators don't denest the roots 
and behave similar as DERIVE’s earlier ver- 
sions: This the screen of the Voyage 200. Tl- 
NspireCAS gives the same results. 


I'FIT^T FT? YTT^T FHT Y F? T Ffir T "l 

L ^—|fllgebra |Ca leather- |Frgm IQ |c lean Up| [ 


■(495 + J495 2 -(1/3) 3 ) +(495 - J495 2 - 

(3 (299 J222 + 4455 )) 1/3 ( "3 [299 J222 > 

3 + 3“ 

>(495 + J495 2 -( 1 / 3 ) 3 ) +(495 - J495 2 - (* 

10.000000 


..<495-JX495 A 2-a/3> A 3>> A <l/3> 


FUNC E/30 


Message 3313: From PAULONE to JERRY GLYNN about VECTOR EXPANSION 

A simple question perhaps about the result of using a vector expression and then simplifying the ex- 
pression. The expression (as authored): VECTOR(n A x , n , -2 , 2) 

When this expression is simplified the result is: 

[2 A x cos (pi x)+2 A x i sin (pi x) , cos (pi x) + i sin (pi x) , ?, l,2 A x] 

I have somewhat of any idea as to what is happening but could you give a simple explanation of what 
Derive is doing with this expression? Positive values of n are okay, of course, but what about the 
negative? 


Message 3314: From JERRY GLYNN to PAULONE about #3313 / VECTOR_EXPANSION 

Let’s simplify this question further so that we can see more clearly what is happening. 

New question: (-2) x If x is 1/2 we have square root of a negative number so Derive simplifies 
(-2) A (l/2) as sqrt(2)*i with the last symbol representing complex number i. Try Simplifying (-2) A x 
and see that Derive is putting this expression into a form that shows the potential complex numbers. It 
is a very weird use of language to say that 2 A x cos(pi*x)+2 A x*i*sin(pi*x) is simpler than (-2) A x. 

Please try out these ideas and come back with a follow on note. Thanks for your question ... this con- 


versation might be helpful to others if made public. 

Recent versions of DERIVE do not show the 
domain of x: 


X 


#1: 

VECT0R(n P n , 

-2 . 2) 

#2: 

x : € Real 



X X 

xl 

#3: 

U-2) , (-1) , 

1, 2 J 

#4: 

x :€ Complex 



I" x-(LN(2) -t jt-l) jt-l-x xl 

*5: 

Le 

, e . 1, 2 J 

#6: 

x :€ Real (-^ p 

0) 


X X 

X xl 

#7: 

U-2) p (-1) ■ 

-■(±1) , 1, 2 J 


May I? 

trig form, but give outputs depending on the 


|Fi7»iriY F£T Y FjT Y FHt Y F£ Y F6t 

| v £— |fl 1 gebra |Ca 1 c |0ther |Prgm 1 0 |C 1 ean 

UF'I ] 

■ seq(n x ? n ? 

-2,2, l) 




tCO£(lT ■ x) 

■ 2 X + sin(n ■ x} ■ 

2 X ■ i 

CO£( JI 

■x) ■+► 

■ seq(n x ? n ? 

-2,2, l) 



2*} 


{( " 1 ) X ’ 2 x ( 

-1) X 

G x 1 

seqCn A x, 

n, “2,2,1 >| 




MAIN 

RAD EXACT 

FUNC 

E/30 



The first calculation was done in Complex 
Format = Rectangular, the second one in 
Complex Format = Real. 

See another difference: DERIVE gives a 
question mark, where the Tl presents 0 X . 


D 
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Message 3333: From KEITH WILLIAMS to PUBLIC about REPLACING COLUMNS IN A 
MATRIX 

I am trying to come up with a user function which will delete a column in a matrix and then put an- 
other column which you specify in it’s place. 

For example a function such as replace_column (m, cl , c2 ) . Where m = matrix and cl = the 
column to be extracted and c2 = the column to be put in the extracted columns place. If anyone can 
come up with such a function I would appreciate it very much. 


Message 3335: From HADUD to KEITH WILLIAMS about REPLACING COLUMNS IN A 
MATRIX 

The following function will replace column no. "nc" with the elements of the vector "c". Is this what 
you want? 


#1: repl ace_col umn (m , nc, c) := VECTOR(VECTOR(IF(jJ = nc, ELEMENT(c, il), 

ELEMENT(m, i I , jJ)), jJ, DIMENSION (ELEMENT (m , 1))), i I , DIMENSION(m)) 


abed 


#2: m := 


e f g h 

i j k 1 . 


#3: repl ace_col umn (m, 3, [p, q, r]) 

a b p d 

#4 : e f q h 

. i j r 1 . 

Now we can write the function in a more compact form: 

repl ace_col umn (m , nc, c) := VECTOR(VECTOR(IF(j_ 


#5: 


= nc, c , m ), j_, 

i_ i_,j_ 


DIM(m )), i_, DIM(m)) 
1 


#6: repl ace_col umn (m, 3, [p, q, r]) 


a b p d 

e f q h 

i j r 1 . 


Message 3338: From KEITH WILLIAMS to HADUD about REPLACING COLUMNS IN A 
MATRIX 

Dear Hadud I am sorry to cause you so much trouble. I had previously left you a message that I could 
not get your solution to my problem to work. It is working now I changed it around a little. This is 
how I changed it. 


BBS : A Lecture via DATA Highway 
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REPLACE_COLUMN(M, NC, C) := VECTOR(VECTOR(IF( j = NC, ELEMENT(APPEND(C) , I), 
ELEMENT(M, i, j)), j, DIMENSION (ELEMENT (M , 1))), i, DIMENSION(M)) 

I really appreciate your help in solving this problem. I tried to figure this out for quite sometime. You 
must be a very smart person. I wish that I could learn to come up with functions as good as you do. Do 
you know of any good books which would teach someone how to better learn to program Derive. 
Again thank you very much for helping me. keith williams 

Message 3340: From JERRY GLYNN to HADUD about REPLACING COLUMNS IN A 
MATRIX 

I agree with Keith's statement that you must be a smart person and I am interested in his question to 
you about how does one learn to write good functions in Derive? Thanks for your help. 

Message 3342: From KEITH WILLIAMS to HADUD about REPLACING COLUMNS IN A 
MATRIX 

Dear Hadud: Your function does work correctly. I just do not understand how the function goes about 
obtaining the results. Could you possibly list the steps that this function goes through and explain how 
the result is obtained. If you could maybe I could better understand how to create user functions. I am 
curious as to how you learned to write them so well. Are you a teacher of math or did you learn what 
you know about writing functions on your own? 

If you could answer these questions I would appreciate it very much. Also any information that you 
could give me to better help me to understand how to write user functions would be greatly appreci- 
ated. keith williams 

Message 3346: From HADUD to KEITH WILLIAMS about REPLACING COLUMNS IN A 
MATRIX 

I'm sorry not to have reacted sooner. I normally log on once a week but I should have realized that 
after leaving a message some questions might come up! 

Looking at the changes you made to make my function work, I suspect the problem was that you were 
in the OPTION-INPUT-(CHARACTER, INSENSITIVE) mode. I assumed the (WORD, SENSITIVE) 
option (my default setting, more about that below) and I should have mentioned that. 

I do not understand why you put in the APPEND function. In the present context it has no effect at all. 
With just one vector-argument APPEND simply returns the argument. Thus APPEND([a,b,c,d]) 
will simplify to [ a , b , c , d ] . Also, in your first message (12/4), the argument following APPEND should 
be "I", not 1, and a final parenthesis is missing. In your second message (12/7) both errors were cor- 
rected. 

Your definition of the function is a good example to demonstrate a sneaky DERIVE quirk as follows: 
Apply your function to the matrix 

abed 

e f g h 

i j k l 

with the 3 rd column to be replaced by [/?, q , r]. 
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This is your result: 



r 

a 

b 

c 

d ' 


REPLACE_COLUMN 


e 

f 

g 

h 

, 3, [p, q, r] 



. i 

j 

k 

1 . 

y 


a b p d 

e f q h 

. 3 2 r 1 . 

The column replacement works as expected but the letters i and j in the matrix are replaced by 3 
and 2. The reason is, of course, that i and j are used as indices within the function. So, when i or j are 
encountered anywhere they are replaced by the current value of the index. Note that this is not an er- 
ror. DERIVE is behaving exactly as advertised. It can be very annoying, however, since the user of a 
function normally does not know which names are used inside the function and therefore should be 
avoided. 

To get around this type of name conflict I always operate in the (WORD, CASE SENSITIVE) mode. 
For function-internal indices and auxiliary functions I use names of at least 2 characters starting with a 
LC character followed by an UC character. This runs contrary to common usage for variable names. 
Thus it is unlikely that a user encounters a name conflict. This explains the use of il and jJ as indices 
in my definition of the function. 

Incidentally, there is no name-conflict problem with argument names used in a function definition, 
since they are treated as dummies. 

As far as gaining some facility with DERIVE programming is concerned, you really don't have to be 
very smart! You might find any introductory text to LISP programming helpful, but the basic princi- 
ples are very simple. The key to LISP (and therefore DERIVE) is the nested function concept. Rather 
than writing sequential statements such as 

y = SIN (x) + A 

z = SQRT (1 + y A 2 ) 

DERIVE wants to see the combined (nested) function SQRT (1 + (SIN(x)+A) A 2) 

To evaluate this expression DERIVE starts with the content of the innermost parentheses (x in this 
case). If a numerical value for x is found it evaluates SIN(x) and then proceeds to the next set of pa- 
rentheses etc. 

Consider now the function REPLACE_COLUMN. The end-product of this function is a matrix. Since in 
DERIVE a matrix is a vector of rows, the outermost (top level) task is to form a vector of rows. So we 
author the expression (let's say it's #1) 

1: VECTOR (row, I , DIMENSION (M)) 
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We haven’t given any thought yet to what ’’row” might be, but we know the number of rows is 
DIMENSION (M), so we let the row index ”1” range from 1 to DIMENSION (M) (the 1 may be omitted, see 
manual). Now we author an expression for ’’row”. We know each row is a vector of elements: 


2 : VECT0R(el em , 3 , DIMENSION (ELEMENT(M , 1) ) ) 


Again we have left open the form of but we prescribe the range of the column index ”J”. Note 

that ELEMENT (M , 1) is the 1st row of ”M”. Its dimension is therefore the number of columns. 

To insert expression #2 into #1 highlight #1 and do a MANAGE-SUBSTITUTE. When the name 
’’row” comes up replace it with ”#2”. This results in the nested expression #3. 

The next step is to define "efcm”. If we substituted ELEMENT(M, I , J) for "gfcm”, the result would be 
the original matrix M, which is OK, except for column no. "AC” which should be replaced. So we 
author 


4: IF(J=NC, ELEMENT(C, I) , ELEMENT(M, I , J)) 


This will simplify to the I-th element of C if we are in the NC-th column. Otherwise we will get the 
original matrix element. 

Now all that's left is to substitute #4 for "efcw” in #3 (giving expression #5) and then name the func- 
tion by authoring 

6: REPLACE_COLUMN(M,NC,C) := #5 

This methodology is called "top-down” (or ”outside-in”) programming. We could equally well have 
used the ”bottom-up” ("inside-out”) approach by recognizing that the innermost task is is the definition 
of an element. Thus we could start with expression #4, embed it in #2 and finally embed the result in 
#1. 

I hope you will find this helpful. After a little practice you will get to the point where you can do all 
these substitutions in your head and write out the entire function in-line (as long as it isn’t too in- 
volved!). 


Final comment from 2008: 


Use the built in function REPLACE(u , v, n) (see Online Help!) 
#7: replace_c(m, nc, c) := REPLACE(c , m\ nc) 1 

a b p d 

#8: replace_c(m, 3, [p, q, r]) = e f q h 

i j r 1 . 


r 

a b 

c 

d ' 



a b 

P 

d ' 


e f 

g 

h 

, 3, [p, q, r] 

~ 

e f 

q 

h 

V 

- 1 j 

k 

1 . 

J 


- 1 j 

r 

1 . 


#9: replace_c 
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Message 3344: From KEITH WILLIAMS to PUBLIC about RECTANGULAR TO POLAR 

Can anyone come up with a user function that will take a matrix such as a 2x1, 3x1 4x1 etc. which 
contains numbers in rectangular form and change them to polar form in degrees and reassemble them 
back as a matrix of whatever size 2x1 etc. that you are working with. 

A function such as Polar(x) := Take like [[4+3i ] , [ 1— 8i ] , [7+3i ]] and change it to polar form 
which would be [ [5 , 36.870] , [8 . 062,-82 . 875] , [7.616,23.199]]. keith williams 

Message 3350: From GREG SMITH to KEITH WILLIAMS about POLAR TO 
RECTANGULAR 

Here are some ideas in response to your messages #3333 and #3344. I have been interested in pro- 
gramming in Derive for some time, especially in connection with implementing procedures in my 
book "Mathematics for Operations Research” which has just come back in print as a Dover Publica- 
tions paperback. I started programming in 1951 when we inherited a Navy computer at George Wash- 
ington University that performed arithmetic using telephone relays, had a magnetic drum, and was 
programmed in octal. It was harder for me to get started in serious programming in Derive than in 
anything I had tried earlier but I must say it has turned out to be supremely rewarding. As a final per- 
sonal note I must note that I am tied up right now in finishing our semester and it will be January be- 
fore I can get back to doing much on this BBS. 

In my experience, the first principle is to start with simple programs and gradually work up to the one 
to be written. The second principle (at least for me because I have concentrated on vectors, matrices, 
and some tensors, as used in my book) is to realize that Derive uses row vectors and not column vec- 
tors (as in my book); this means that even I, who for decades have forced students to use one-rowed 
matrices to represent, for example, Frechet derivatives, must pay due respect to row vectors. Maybe 
the next principle is to check every individual command (or program) to see that it works correctly. 

With regard to #3333, consider starting as follows to do the job in three steps, each a small step to- 
wards the goal. 

REPLACE_VECTOR_ELEMENT (v , j 0 , new) : =VECTOR(IF(n_=jO , new, ELEMENT(v , n_)) , n_, 

DIMENSION(v)) 

REP LAC E_M ATRIX_R0W ( a , i 0 , v) : =REPLACE_VECTOR_ELEMENT (a , i 0 , v) 
REPLACE_MATRIX_COLUMN_BY_VECTOR(a , i 0 , v) : =REPLACE_MATRIX_ROW(a 1 ,i0,v) 1 

You may wish to use abbreviated names but note that "BY VECTOR” in the preceding name is in- 
tended to stress that we are using a (Derive) vector and not a one-row or one-column matrix. Try these 
to see that they work. 

With regard to #3344, the following should indicate how you could move toward this goal. First, con- 
sider how you could replace one element in a matrix as the counterpart of the first command used 
above. In order to keep it simple, it might help to use the following two commands which are cer- 
tainly not necessary but they can aid the eye in working with matrices; I would not use these for (De- 
rive) vectors that were not actually matrices. 
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WIDTH (a) : =DIMENSION(ELEMENT(a , 1) ) 

DEPTH(a) :=DIMENSION(a) 

REPLACE_MATRIX_ELEMENT (a , i 0 , jO , new) : =VECTOR(IF(m_=i 0 , VECTOR(IF(n_=j 0 , new, 

ELEMENT(a, m_, n_)) , n_, WIDTH (a)) , ELEMENT(a, m_)) ,m_, DEPTH(a)) 

Now, finally, to address your #3344, try to adapt the following to your needs. The simplest command 
to replace one matrix by another whose elements are functions of those for a given matrix might be the 
next one. 

COPY_MATRIX(a) :=VECTOR(VECTOR(ELEMENT(a,m_, n_) , n_, WIDTH(a)) , m_, DEPTH(a)) 

Check this to see that it works. Now, simply Replace 'ELEMENT (a , m_, n_)' in the former command 
by 'V(ELEMENT(a, m_, n_))' and assign the name 'SQRT_ELEMENT(a)'. 

SQRT_ELEMENT (a) : =VECTOR(VECTOR(SQRT(ELEMENT (a , m_ , n_) ) , n_ , WIDTH (a) ) , m_ , DEPTH (a) ) 

Check this to see that it works on matrices with nonnegative elements and then move on to your case 
of interest. Good luck, w.h.marlow 

Message 3354: From KEITH WILLIAMS to W.H.MARLOW about RECTANGULAR TO 
POLAR 

Dear MR. MARLOW I was looking over your functions which you left in a file for me. I came up 
with a function ALL (M , R1 , R2 , R3) : = where M= MATRIX, 3x1 . R1 to R3 are the rows of the 3x1 ma- 
trix which will be changed into polar form. I uploaded a file named POLAR.ZIP in the file area so that 
you could download it. This function will just work for a 3x1 matrix. I haven't figured out how to 
come up with a function which will take a matrix of 2x1, 3x1, etc.... and change the rows from rectan- 
gular to polar form. Maybe you can look at the file that I uploaded come up with a function which will 
do this. Without your help and the help of HADUD I could never have come up with this function. So 
I guess I am getting better at figuring out functions. Look at the file which I uploaded and maybe you 
can come up with the function POLAR(M):= which will take any size matrix 2x1 and up into polar 
coordinates. Thank you very much for your interest in my problem. KEITH WILLIAMS 

Message 3355: From HADUD to K WILLIAMS about RECTANGULAR TO POLAR 
CONVERSION 

Allright, let’s try it once more! This time we’ll use the bottom-up approach. 

What you want to do is to take the first (and only) element of the I-th row of your (n x l)-matrix and 
convert it to polar form. Then you want to build up a result matrix from this. 

OK, you get the element by authoring 


1: ELEMENT(M, I ,1) 






D 


-N-L#17 


BBS: 


A Lecture via DATA Highway 


P 13 


To convert any complex number Z you have to compute 2 numbers: amplitude and phase. To do that 
in one swoop, you make a vector out of the results by authoring 
2: [ABS(Z) ,180/pi PHASE(Z)] 

Now highlight #2, do a MANAGE-SUBSTITUTE and when Z comes up replace it with #1. The result 
is 

3: [ABS(ELEMENT(M, 1,1)) ,180/pi PHASE(ELEMENT(M , I , 1)) ] 

Now you have an expression for the I-th row of your desired output matrix. To put all the rows in a 
matrix and create the desired function you author 

4 : REC2 POL (M) : =VECT0R(#3 , 1 , DIMENSION (M) ) and that’s it! 

Note that I used single-character names throughout, so this will work in the OPTION-INPUT- 
(CHARACTER, INSENSITIVE) mode. To bring the result closer to customary notation, you might 
want to replace the comma in expression 2 with a <, which can be read as an angle sign: 

2: [ABS(Z)<180/pi PHASE(Z)] 

DERIVE now interprets what’s inside the brackets as a relational expression, but that's OK as long as 
don't do any further processing of the result matrix. 

Message 3356: From KEITH WILLIAMS to HADUD about RECTANGULAR TO POLAR 

Dear Hadud You have come to my rescue one more time. What can I say to you except that I really 
do admire your understanding of mathematics. You must be a genius. I sure wish I could figure out 
problems like you do. I bet that there are not many problems in math that you cannot figure out. If I 
need anymore help with any problems I know who to consult. The reason that I was interested in ma- 
nipulating matrices is because I am taking electrical courses that involve using matrices to solve elec- 
trical circuits. In ac circuit analysis you have to deal with rectangular and polar coordinates and phase 
angles. I really appreciate your help on my problems. Maybe I can help you sometime. Thank you 
very much for your help Hadud. keith williams 

Message 3357: From KEITH WILLIAMS to HADUD about POLAR TO RECTANGULAR 
CONVERSION 

DEAR HADUD, The function that you gave me to change from rectangular to polar helped me out a 
lot. I took your function and also made it where it would change from polar to rectangular. 

pol 2 rec(m) :=vector([element(m, i , l)*cos (element (m, i , 2 )* (pi / 180 ))+ 

element(m,i , 1 ) * (si n (element (m, i , 2 )* (pi / 180 ))*#i)] ,i ,dimension(m)). 

This will change from polar to rectangular, keith williams 

Message 3358:From HADUD to KEITH WILLIAMS about POLAR TO RECTANGULAR 
CONVERSION 

Thanks for communicating your neat P0L2REC function. I get the impression that you are now well on 
your way to becoming an expert in DERIVE functionality! 

If you prefer to work in the character input mode I have one suggestion: For the reasons outlined in 

my message #3346, use letters followed by underscore(s) for indices inside your function (I_ or I in 

place of I in the function P0L2REC). DERIVE allows this type of variable names even when you are in 
the character mode (see version 2 manual , p.38) 
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Message 3359: From KEITH WILLIAMS to HADUD about #3358 / POLAR TO 
RECTANGULAR CONVERSION 

Dear Hadud I know that you have heard enough about matrices and polar to rectangular conversion, 
but I am going to ask you one more problem. You came up with the function rec2pol which would 
change a 2x1, 3x1 etc ... matrix from rectangular to polar form. I got to thinking that I would try to 
come up with a function which would take a 2x2, 3x3, 4x4 etc... matrices which contained elements in 
rectangular form to polar form. Well, I am still working on it. I tried to use your reasoning on how to 
figure it out using the methods which you used in figuring out the other problems. So far I haven't had 
much luck. I know that the end product of this function is a matrix. So think that you would need vec- 
tor(row, i , dimension(m)). I got this from one of the problems that you solved. Then vec- 
tor (el em, j , dimension (el ement(m,l))). And then replace row with the last expression, which 
would yield vector (vector (el em, j , dimension(element(m,l))) ,i ,dimension(m)). I know 
that something needs to be put in place of elem, but I have not figured it out yet. I tried putting 
[ abs (el ement (m , i , 1) ) , 180/ pi phase (el ement (m , i , 1) ) ] but this replacement doesn't seem to 
work, can you help me with this problem. If you will I will try my best to figure out these problems 
concerning matrices myself. 

Thank you keith williams 

Message 3360: From HADUD to KEITH WILLIAMS about RECTANGULAR TO POLAR 
CONVERSION 

Everything you did is perfectly OK. The only problem is that, when replacing "e/em", you forgot that 
you are no longer dealing with a (n x 1) but with a (n x m) matrix. Your element references therefore 
should be "el ement (m , i , j)", not "el ement (m , i , 1)". If you make that change you will get a work- 
ing function but the display will not be what you want. The reason is as follows: 

Remember that DERIVE works inside-out. Thus when it encounters the innermost brackets [] it inter- 
prets their contents as a row. The next VECTOR function (which has the index J) then forms a (m x 2) 
matrix. The outermost VECTOR function (index I) then forms a vector of matrices, i.e. it will display n 
(m x 2)-matrices side-by-side. 

The net result is something that looks like the transposed of the desired matrix. This is another 
DERIVE quirk (not an error) that can be frustrating. 

There are two ways around the problem: 

1) Put another set of brackets around the replacement for "efcm": 

[ [ABS(ELEMENT(M ,1,3)), 180/pi PHASE(ELEMENT(M ,1,3))]] 

Now " elem " is a (1 x 2) matrix and the final result will be a matrix of matrices with the desired ap- 
pearance. 

2) Omit the brackets and replace the comma with a relational symbol: 

ABS(ELEMENT(M, I , J)) < 180/pi PHASE(ELEMENT(M , I , J)) 

Now "elem" is just an expression, not a vector, and the result will be a (n x m) matrix of expres- 


sions. 
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Message 3361: From KEITH WILLIAMS to HADUD about RECT.TO POLAR CONVERSION 

Hello Hadud Your function worked perfect again. At one time I had the function down just like you 
had it except for the extra brackets. While we are really into matrix manipulation I have one more 
problem for you since you have not failed yet to come up with a solution to every problem which I 
have burdened you with. May I ask you one more question. I would say that this would be the last one 
but I would not want to lie to you. 

You really amaze me how you can figure out all of these matrix problems. What I would like to ask 
you is that your last function would convert a 2x2, 3x3 etc., matrix from rectangular to polar coordi- 
nates but can you develop a function which would do the opposite: convert a 2x2,3x3,ect... from polar 
form to rectangular form? If you can I don't think that I have any more questions about matrix manipu- 
lation at the present time. Thank you very much for all the help which you have given me on this sub- 
ject. keith williams 

Message 3363: From KEITH WILLIAMS to HADUD about MATRIX REDIMENSION 

Dear Hadud I have a hp48 calculator and I was looking through the manual which came with it on the 
functions which would manipulate matrices. In the manual there was a function called redimension 
which would take a 2x2, 3x3 ect.. and redimension it. For example in the book there was an example 
which showed a 3x3 matrix with the element [[1,2, 3], [4, 5, 6], [7, 8, 9]] and was redimensioned to a 
9_row by l_col matrix with the elements redimensioned as [[1], [2], [3], [4], [5], [6], [7], [8], [9]]. Would a 
function such as this be very hard to develop for derive? A function such as rdm(m, row, col) : = 
where m = matrix row = number of rows and col - number of columns. 

I was just curious wheather or not a function like this could be developed for DERIVE. Thank you for 
taking an intrest in all these problems which I have asked you concerning matrix manipulation, 
keith williams 

Message 3364: From HADUD to KEITH WILLIAMS about RECT.TO POLAR CONVERSION 

I’m a little disappointed that you couldn't do the POL2REC function yourself after all my explana- 
tions! Here it is: 

P0L2REC(M) : =VECTOR(VECTOR(ELEMENT(ELEMENT (M ,1,3) ,1,1)* (COS(ELEMENT 

(ELEMENT(M , I , J) , 1 , 2) *pi /180)+#i *SIN(ELEMENT(ELEMENT(M , I , J) , 1 , 2) *pi /180)) , 

3 , DIMENSION (ELEMENT(M , 1))) , I , 1 , DIMENSION(M)) 

This assumes the double-bracket version of REC2POL, i.e. 

REC2P0L(M) :=VECT0R( VECTOR ([ [ABS(ELEMENT(M ,1,3)), 180/pi *PHASE(ELEMENT(M ,1,3))]], 

3 , DIMENSION (ELEMENT(M , 1))) , I , DIMENSION (M)) 

The redimension function is best done recursively by introducing the auxiliary function APPV_: 

APPV_(n ,M) :=IF(n=l, ELEMENT(M,1) , APPEND (APPV_(n-l , M) , ELEMENT(M , n))) 

REDIM(M) :=ELEMENT([VAPP_:=APPV_(DIMENSION(M) ,M) , VECTOR([ ELEMENT(VAPP_, I) ] , 

I , DIMENSION (M)*DIMENSION( ELEMENT(M,1))) ] , 2) 


All of the above functions will work for any matrix M, square or rectangular. 
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Message 3366: From KEITH WILLIAMS to HADUD about RECT.TO POLAR CONVERSION 

Dear Hadud I did come up with a function the way I could understand it. Mine is I think like yours 
except that it is longer. It is as follows: 

pol2rec(m) : =vecto r (vector (element (element (el ement(m,i , j) *1,1) *1,1)* 

cos (element (el ement(m,i , j) ,1,2)* (pi /180) )+element (element (el ement(m,i , j)*l,l) 
*l,l)*sin(element(element(m,i , j) ,l,2)*(pi /180))*#i , j , dimension 
(element(m,l))) ,Zi ,dimension(m)) 

keith williams 

Message 3367: From KEITH WILLIAMS to HADUD about MATRIX REDIMENSION 

Hadud the functions which you gave to me the other day worked. I guess I didn't make myself clear as 
to the kind of function I wanted. The example in the book took for example a 3x3 matrix and redimen- 
sioned into like a 3x1, 1x3, 1x9, 9x1, 2x4, 4x2, etc., just whatever dimension that you could assemble 
out of what ever size matrix which you were working with. The example took a 
[[1,2, 3], [4, 5, 6], [7, 8, 9]] matrix and changed it into different sizes, a 9_row l_col would be 
[[1], [2], [3], [4], [5], [6], [7], [8], [9]] A 3x2 would be [[1,2], [3, 4], [5, 6]] A 2x3 would be [[1,2, 3], [4, 5, 6]]. 

I have been trying to come up with a function which would do this. 

For example redi m_mat ri x(m , col , row) : = where m = matrix ,col = number of columns wanted, and 
row=number of rows wanted. I have put to together some of it but it will just take the first element of 
the matrix and make whatever size matrix which you specify with all the elements in the new matrix 
of element 1 .What I have is 

Redi m_matrix(m, col , row) :=vector (vector (el ement(m, 1,1) , j ,l,col) , j ,1, row) . 

This will recreate any size matrix which you can redimension out of the original matrix but every ele- 
ment win the matrix will be the same as element number one. I just thought that I would show you 
what I have come up with so far. Look at what I have come up with and tell me if I a proceeding in the 
right direction. If you come up with a function which will do this let me know what it is. Mean which I 
will keep trying to figure this one out. Ill be waiting to hear from you. Maybe I can figure this one out 
and show you that I have at least learned something. Thank you keith williams 

Message 3368: From KEITH WILLIAMS to HADUD about MATRIX REDIMENSION 

Dear Hadud I just wanted to show you this function which I came up with. It is not anything like what 
you can do but at least its a beginning. This function called 

Matrix_reduce(m,dw,dd) : =vector (vector (el ement(m,i , j) , j , dimension (element (m,l))- 
dw) ,i ,dimension(m)-dd) . 

This function will take any size matrix reduce the number of rows or columns by the amount which 
you specify, dw = deduct width, dd = deduct depth, m = matrix. I just wanted to show you this, keith 
williams 
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Message 3369:From KEITH WILLIAMS to HADUD about #3367 / MATRIX REDIMENSION 

Dear Hadud I am still trying to come up with a function which will redimension a matrix into different 
dimensions. I made a little progress today I came up with 

Matrix_redim(m,col , row) : =vector (vector (el em, j ,col) , I, row) . 

This will create any size matrix which you specify by col and row with the word elem in place of the 
elements of the matrix. I haven't figured out what to put in place of elem to make it redimension the 
elements of the matrix that you would be working with. I'll keep trying I know without a doubt that 
you can figure it out. If you do let me know, keith williams 

Message 3371:From KEITH WILLIAMS to HADUD about RECT.TO POLAR CONVERSION 

Dear Hadud The function which you gave me to redimension a matrix into a l col by x rows works. I 
figured how to make the function work faster. 

The one that you gave me took a 20x20 matrix and changed it into a l col by 400_row matrix in 
309_sec. The way which I came up with did in in .3_sec. The function is: 

redim2(m) :=extract_l_row(adjoi n_element (append (m) , m) ,1). 

I made this function from the following: 

ext ract_l_el ements (m , e) : =[ el ement (m , e) ] 

Ext ract_l_row(m , cl) : =ext ract_l_el ements (m , cl) ' 

Adjoi n_element(ml,m2) :=append([ml] ,m2) keith williams 

Message 3372: From HADUD to KEITH WILLIAMS about REDIMENSIONING A MATRIX 

Bravo! Your REDIM2 function is indeed much more elegant than my monstrosity! I guess I can retire 
as your teacher! I learned from your function that you can use APPEND (M) to convert a general matrix 
M to a single vector, which is something I never realized. Given that, I can still do you one better, 
though! Try 

REDIMl(M) :=[APPEND(M) ] ' 

This also holds the key to your general REDIM function: 

REDIM(M , r , c) : =VECTOR(VECTOR(ELEMENT (APPEND(M) , (I-l)c+J) , 3 , c) , I , r) 

This will generate an (r x c) matrix from the elements of the arbitrary matrix M. I hope I've understood 
the problem correctly now! 

Message 3373: From KEITH WILLIAMS to HADUD about REDIMENSIONING A MATRIX 

Dear Hadud I'm the one who should say bravo! to you. The redi m(m , r , c) function was just what I 
wanted. I don't understand it but it works exactly as it should. Your rediml function was faster than 
my redim2 function which just goes to show you that you are still a lot better than I am at figuring out 
functions. I can't think of any more questions concerning matrices, but don't retire yet as my teacher. I 
want to say to you that I appreciate all your help on all these problems which I as sure were not any 
challenges to you. I have enjoyed communicating with you over this bulletin board during the past few 
weeks. Again I want to thank you for taking an interest in my problems in math. Your friend and ad- 
mirer keith williams 
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Solving the problems with recent version of DERIVE: 

redim(m_, r_, c_) VECTOR(VECTOR((APPEND(m_)) 

# 1 : 


, k , O , i , O 


(i - l)-c_ + k 

redim 2 (m_, r_, c_) := 

If r_* c_ t DIM(APPEND(m_)) 

# 2 : "wrong dimension" 

VECTOR(VECTOR( (APPEND (m_)) 4 , ((i - l)*c_ + k), k, 1, c_) , i, 1, r_) 
#3: m := VECTOR([ RANDOM(IOO) , RANDOM(IOO) , RANDOM(IOO) ] , k, 4) 


#4 : m 


' 33 

58 

23 ' 

51 

93 

96 

57 

91 

90 

. 54 

83 

44 . 


#5: redim(m, 2, 6 ) 


33 

58 

23 

51 

93 

96 ' 

57 

91 

90 

54 

83 

44 . 


# 6 : redim(m, 6 , 2 ) 


33 58 

23 51 

93 96 

57 91 

90 54 


83 44 

33 58 23 

51 93 96 

57 91 90 

# 8 : redim2(m, 3,3)= wrong dimension 


#7: redim(m, 3, 3) 


# 1 : ml := 


4 + 3 • l 

1 - 8* L 

7 + 3 • L 


#2: REC2P0L(M) VECTOR^VECTOR^ |ELEMENT(M, i, j)|, 


180 


■ PHASE (ELEMENT (M , 


i > j)) 


, j, DIMENSION (ELEMENT(M , 1)) 


), i, 


DIMENSION(M) 


#3: REC2P0L(ml) 

simplified and then approximated: 

[5, 36.86989764] 
[8.062257748, -82.87498365] 
L [7. 615773105, 23.19859051] 


# 5 : 
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#8: P0L2REC(m2) 

simplified first and then approximated: 

' 4 + 3-l 

#10: 1 - 8-l 

. 7 + 3-l 

Version for DERIVE 6 

180 

|v|, (PHASE(v)) |, v, m_ 

TT 1 

#12: REC2P0L_(ml) 

Expression #12 approximated 

5 36.86989764 

#13: 8.062257748 -82.87498365 

. 7.615773105 23.19859051 

5 36.86989764 

#14: m2_ := 8.062257748 -82.87498365 

. 7.615773105 23.19859051 

P0L2REC_(m_) := VECT0R(v_ -(C0S(v_ *1°) + l-SIN(v_ *1°)), v_, m_) 
#15 : 12 2 

#16: P0L2REC_(m2_) 

#16 approximated: 

#18: [4 + 3-l, 1 - 8-l, 7 + 3-l] 


#11: REC2P0L_(m_) := VECTOR 


More Forum contributions on page 50 
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Ebene Algebraische und 
Transzendente Kurven (6) 

Thomas Weth, Wurzburg, Germany 

Cassinische Kurven - Die Lemniskate von Bernoulli 

Im Jahre 1680 beschrieb Johann Dominicus Cassini eine Kurvenschar, die bis heute mit dem Namen ihres Er- 
finders verbunden blieb. Uber den Anlass ihrer Entdeckung schreibt sein Sohn Jakob Cassini im Jahre 1749 
sinngemaB: „Durch genaue Beobachtung der Sonne fand mein Vater eine von der Ellipse verschiedene Kurve, 
die dazu dient, die exakte Bewegung der Sonne und ihre verschiedenen Entfemungen von der Erde zu beschrei- 
ben.“ Obwohl sich herausstellte, dass die Cassinischen Kurven in keinerlei Beziehung zum eigentlichen 
Ausgangsproblem stehen, dienten sie doch immer wieder als Gegenstand mathematischer Untersuchungen - 
insbesondere um die Tragfahigkeit algebraischer Methoden und der Infmitesimalrechung auszuloten. 
Aufmerksamkeit schenkte man einem Spezialfall, der Bemoullischen Lemniskate, mit deren Rektifikation (Lan- 
genmessung) sich Graf Fagnano eingehend beschaftigte. In geeigneter Weise entwickelt und verallgemeinert, 
bilden seine Uberlegungen den Kern der Theorie der elliptischen Funktionen. 

1680 J.D. Cassini decribed a family of curves, connected with his name until today. His son, Jakob 
Cassini, wrote in 1749: “...by exact observation of the sun may father found a curve, different from an 
ellipse, which can serve to describe the exact motion of the sun and her different distances from 
earth." Although it was found out that the curves of Cassini don't have any connection with the starting 
problem, they have been used as an object of mathematic investigation, especially to check the valid- 
ity of algebraic methods of calculus. Special attention was paid to a special case: the Lemniskate of 
Bernoulli. Count Fagnano dealt with its rectification and his ideas can be seen as the foundation of the 
theory of elliptic functions. 

Definition und Konstruktion der Kurve 

Gegeben sind zwei Punkte Bi und B 2 (Brennpunkte) mit der Entfernung 2c. Die Punkte P, bei denen 
das Produkt ihrer Entfemungen pi und p 2 zu den Brennpunkten einen festen Wert a 2 besitzt, bilden die 
Cassinischen Kurven: pi • p 2 = PBi • PB 2 = a 2 . 

Given are two points Bi and B 2 (focal points) with = 2c. The curves of Cassini are the 
sets of points P, whose product of the distances to the focal points pi • p 2 = PBi • PB 2 = a 2 
(const.). See fig. 3. 

Legt man die Brennpunkte symmetrisch zum Urspmng eines kartesischen Koordinatensystems auf die 
x-Achse (vgl. Abb. 3), so erhalt man in kartesischen Koordinaten: 

■\j ' y 2 + (x + cf y jy 2 + (x - cf = a 2 und durch Quadrieren: (x 2 + y 2 f - 2 c 2 (x 2 - y 2 ) = a 4 - c 4 . 

Zur Gleichung in Polarkoordinaten kommt man durch Einsetzen von x = r cos ^undj; = r sin cp (mit 
DERIVE): r 4 - 2c 2 r 2 cos 2cp = a A - c 4 . Lost man diese Gleichung nach r auf und plottet in Polarko- 
ordinaten, so liefert DERIVE: (soLve for r and plot in polar co-ordinates:) 



Abbl: Cassinische Kurven mit DERIVE 
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Mit dem Schieberegler lasst sich nun auch die implizite Darstellung treffen. Daneben ist das Bild einer 
Schar von Cassinischen Kurven in impliziter Form gegeben. 


Applying slider bars for a and c and the implicit form it is easy to investigate the curves. On 
the right hand side you can see a family of "Cassinis". 






Abb2: Cassinische Kurven mit Schieberegler und impliziter Plot (Kurvenschar) 

Cassini Curves with slider bars for a and c and a family of implicit plots 


Zur Zirkel- und Linealkonstruktion der Kurve 
benotigt man zu gegebenem a 2 zwei Strecken pi 
und p 2 , fur deren Produkt gilt: pi • p 2 = a 2 . Um 
die gegebenen Brennpunkte Bi und B 2 werden 
dann Kreise mit diesen Radien gezeichnet. Die 
Schnittpunkte dieser Kreise sind dann Punkte 
der Cassinischen Kurven. 

Die Struktur der Gleichung pi • p 2 = a 2 legt zur 
Konstruktion von pi und p 2 die Verwendung 
des Kathetensatzes (oder des strukturgleichen 
Hohensatzes) im rechtwinkligen Dreieck nahe. 

Gibt man auf der x-Achse den Punkt A(a,0) vor, 
und verbindet A mit einem beliebigen Punkt R 
der y-Achse, so ist das Dreieck AOR rechtwinklig. Die Hohe von O auf die Hypotenuse AR liefert den 
HohenfuBpunkt T und den zu Kathete OA zugehorigen Hypotenusenabschnitt AT der Lange pi. Mit 
der Hypotenusenlange p 2 gilt nach dem Kathetensatz: pi • p 2 = a 2 . 

Zeichnet man nun die Kreise k(Bi, pi) und k(B 2 , p 2 ), so bilden deren Schnittpunkte P - soweit sie exis- 
tieren - Punkte der Cassinischen Kurven (siehe Abb. 3). 

To construct the curve using ruler and compass you need to a given a 2 two distances pi and 
p 2 with pi • p 2 = a 2 . Then you have to draw circles with centres B 1j2 and pi j2 as radius. The 
intersections of these circles are points of the curves. The structure of pi • p 2 = a 2 leads to 
use the Pythagorean Theorem. Take any arbitrary point A(a,0) on the x-axis and connect A 
with an arbitrary point R on the y-axis, then the triangle AOR is a right one. The height OT 
with respect to AR gives T, which divides the hypotenuse AR = p 2 . With AT = pi you find fi- 
nally RA • TA = OA 2 = pi • p 2 = a 2 . So draw two circles k(B l5 p^ and k(B 2 , p 2 ), and their inter- 
section points are - if they do exist - two points of the curves. 

See figure 3. 



Abb. 3: Zur Konstmktion der Cassinischen Kurve 
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Die Bernoullische Lemniskate 

Der Sonderfall c = a liefert einen Spezialfall der Cassinischen Kurven, die Lemniskate von Bernoulli, 
deren Namen sich ableitet vom griechischen A,8 |uvigko(^ (Schleife, Band). 


The special case c = a gives the Lemniscate of Bernoulli. (^spviaKo^ = loop, cord) 

In bipolaren Koordinaten ergibt sich pi • p 2 = c 2 (bipolar co-ordinates). 

Als algebraische Kurvengleichung erhalt man daraus: (x 2 +^ 2 ) 2 = 2c 2 (x 2 -y 2 ) 
und in Polarkoordinaten: r 2 = 2 c 2 cos 2cp. 



I2AF ■ I"COP 1 

J 

+ ^ 

,'rv 

J 

1 1=1. 40cm 
\2=0.53 

. y . . * 

\ 

\, 

O.S6cm 

0.86cm 
a E = 0. 74 

MAIN DEG EXACT FUNC 


Abb. 4. Lemniskate von Bernoulli (with Cabri on the PC and on the TI-Handheld) 

Eine weitere Konstruktion, die sich als „wortliche Ubersetzung“ der algebraischen Kurvengleichung 
verstehen lasst findet man in „Kurven in geometrischer und algebraischer Sicht“ (Didaktik der Ma- 
thematik 2, 1992 (112-138)). 

Einige interessante Eigenschaften der Lemniskate, die sich mit DERIVE iiberprufen lassen: 


• Die Kurve schneidet sich im Ursprung selbst in einem rechten Winkel. 

• Die Lemniskate ist eine anallagmatische Kurve, dh sie wird unter einer (geeigneten) Inversion am 
Kreis auf sich selbst abgebildet. 

• Den von der Lemniskate eingeschlossenen Flacheninhalt erhalt man (DERIVE !): 

71 

A = 4(- j r 2 dty ) = 2c 1 
2 <p=o 

Der eingeschlossene Flacheninhalt ist also halb so groB wie ein Quadrat, dessen Seitenlange gleich der 
Entfemung der beiden Brennpunkte ist. 

The lemniscate shows some interesting properties, which easily can be proofed using 
DERIVE: 


• The curve is intersecting itself orthogonally. 

• The lemniscate is anallagmatic, i.e it can be mapped by an appropriate inversion at a cir- 
cle onto itself. 

• The area enclosed by the lemniscate is c 2 , is half the area of a square with the distance 
between the foci as side length. 


D-N-L#17 


Thomas Weth: A Lexicon of Curves (6) 


P 23 


Details from the DERIVE session: 


#16: 

4 2 

r - 2-c • 

#17: 

4 

S0LVE(r - 

#18: 

r = - 7^2 - | . 

#19: 

r = V2-|c| 

Area: 

0 < (p < 7t/ 4 g 

#20: 

r < J2 - | c | 


2 2 



#22: r < J2- |c| -7CC0S (2 .()>)) 


#23: VECTOR 


U2- |c|.y(C0S(2.<t>)), <p], <p, 0, , 

4 5000 / 



#21: 4- 


n/4 


L.J 

2 0 


(V2.|c|.7(COS(2.<p))) d<p 


= 2 • c 



Family Cassini (Memories on DOS-Times) 
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A DERIVE-Toolbox 

Eugenio Roanes-Lozano , Eugenio Roanes-Macias & Josef Bohm 

This contribution has a nice history. Eugenio jr sent a FAX with his (( A Trick about 
2D Plots in DERIVE (( and I gave the answer, that I had written a little utility to draw 
broken lines and shaded areas to complete Sergey 's labelling utility from DNL#16. So 
Eugenio suggested publishing a Spanish - Austrian co production. This is the result: 

A Trick About 2D Plots in DERIVE 

Abstract 

There are some commands in DERIVE that have to be selected from the Command line and do not 
have a correspondent function in the Author line. This is the case with colours. It is not possible to 
choose directly the desired colour for the plot of a function from the Author line. 

This is sometimes an inconvenience. For instance, when preparing a .MTH - file for a Teaching 
Unit, if we want to identify a certain curve by its colour. 

We shall show how it is possible to choose the colour of a plot from the Author line (and to keep 
it unaltered even if it is plotted more than one time). It works both in DERIVE 2 and the new DERIVE 
3. 


1.1 Plotting a function 

We shall assume that our hardware uses a 16-colour VGA graphics-card. If that is not the case, what 
follows also stands, altering it accordingly to the number of colours of the corresponding graphics- 
card. 


In DERIVE, by default, if the last plot was drawn in colour /, the next plot will be drawn in col- 
our i-1. The other possibilities given in the Command line are: 

• From now onwards, plot all the functions in the same colour 

(Plot/Options/Color/Auto/No) 


• Choose manually the colour of the next plot 

(Plot/Options/Color/Plot/Next Plot/ ...) 

If we leave it as it is by default, after 15 plots (the background colour is not used), the colour of the 
plot will be the same. Moreover, a vector can be plotted. Therefore we can merge these two ideas to 
obtain any plot in the desired colour. We have to build a function DRAW(fu,CO) that constructs a 15 
items vector, with the function fu in position 15 - CO+1 . The other 14 functions must plot nothing (in 
the not required colours). We have chosen y = T, but this can obviously be changed. 


Comment on this for the recent version of DERIVE: 
VGA-graphics adapters are the past. The order of 
the colours did also change. We can choose from 16 
colours as ever. The standard colours are presented 
in a box and can be numbered from 1 to 16 (1 = red, 
2 = blue, ..., 16 = white): 


Display Options 


Axes | Cross | Grids | Points Color | 
Plot 


Next color: 


Background 


!l^ Automatically chanc|B|® ® slots 

mrur 

mmmr 
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Let's define the function DRAW(fu,CO) that draws the function fu in colour CO this way: 

#1: draw(f, co) := APPEND(VECTOR(l , k, co - 1) , [f], VECTOR(l, k, 15 - co)) 

Instead of ’’plotting” the imaginary unit you can plot ’’nothing” expressed by the questions mark! 

#2: draw(f, co) := APPEND(VECTOR(? , k, co - 1) , [f], VECTOR(?, k, 15 - co)) 

In order to have the colours in a correct ordering, before using DRAW for the first time, it is conven- 
ient to start a new DERIVE session to confirm in the Display Options Color Box, that red is the colour 
for the next plot and that "Automatically change color of new plots" ticked. 

Let's create the first ’’programmed” graphs of functions: 


#3: Ldraw(x , 5), draw(- x , 4) , draw(SIN(x), 1), draw(C0S(2-x) , 3)J 

Expression #3 should show the first parabola in purple (colour #5), the second one in light green (col- 
our #4), the sine-wave in bright red (colour #1) and the cosine in black (colour #3). After each plot the 
colour table is set with red to start. 



You see that this is working properly. 


1.2 Creating a Colour Table 

Assuming that your colour palette looks like the palette shown on page 24 we can assign: 

#4: [red := 1, blue := 2, black := 3, l_green := 4, purple := 5, m_green := 6, 

d_gray := 7, l_blue := 8] 

#5:= [grayblue := 9, l_gray := 10, d_blue := 11, yellow := 12, green := 13, 
d_brown := 14, pink := 15, white : = 16] 




1 




# 6 : 


draw(TAN(x) , l_bl ue) 


p 26 


A Spanish-Austrian DERIVE - Toolbox 


D-N-L#17 


1.3 Plotting several functions at the same time 

Now, to plot several functions with just one line of code, we shall do a Plot of a VECTOR of the vec- 
tors created by DRAW (i.e., a Plot of a matrix). 

2 

#7: VECTOR(a-x , a, 4) 

This does not address the colours. Assume that you want to have all parabolas plotted dark brown, 
then enter and plot: 


2 

#8: VECTOR(draw(a-x , d_brown) , a, 4) 



Plot in different colours y with j ranging from 1 to 4. 

j 

#9: VECTOR(DRAW (x , j),j,l,4) 



Final example: draw a green triangle 

#10: [a := [1, 1], b := [6, 1], c := [3, 4]] 

#11: draw([a, b, c, a], green) 



References: 


[1] J.L.Llorens: Introduction al uso de DERIVE , Universidad Politecnica de Valencia, 1993 

[2] , [3] A.Rich, J.Rich, D.Stoutemyer: DERIVE User's Manual 2 & 3, Soft Warehouse 
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In the last DNL you could find Sergey 's Labelling Tool. Often it is useful to draw bro- 
ken or dotted lines and to shade areas in a plot to demonstrate or to emphasize a cer- 
tain fact. Using the next little functions you can easily do this and draw pure DERIVE 
- sketches without using any other drawing program. 

Broken lines and shaded areas 

See first four applications of my functions: 







Figure 4a (DERIVE 6 - shading using inequalities, plotting u < y < v) 
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I used the following DERIVE Utility file BROKEN.MTH: 
#1: [Notation Decimal, Notati onDi gi ts := 3] 


NORM(xl, yl, x2, y2) := [x2 - xl, y2 - yl] 

#2: 2 2 

7((x2 - xl) + Cy2 - yl) ) 


#3: STR(xl , yl, x2, y2, 1, p) := APPEND VECTOR| [[xl, yl] + (i - 1) • (1 + p) 


•NORM(xl, yl, x2, y2) , [xl, yl] + (i-1 + i-p - p)-NORM(xl, yl, x2, y2)], 


( 2 2 ' 

VCCy2 - yl) + (x2 - xl) ) 

i, 1, FLOOR , [xl, yl] + 

>. p + i JJ L . 


2 2 2 
/(xl - 2*xl*x2 + x2 + (yl - y2) ) 

1 + P 


(1 + p)-NORM(xl, yl, x2, 


y2), [x2 , y2] 


#4: ORDNER(x, y, 1, p) := [STR(0, y, x, y, 1, p) , STR(x, 0, x, y, 1, p)] 


#5: SHADECu, v, x, a, b, n) := VECTOR 


i • (b - a) 


i • (b - a) 


1 im v 

x->a + i • (b - a)/n 


1 im 

x->a + i • (b - a)/n 


, i , 0 , n 


#6: SHADE_C(u, v, x, a, b, n, ca, ce) APPEND([ [IF(ca < x < ce, u, ?)]], 

[[IF(ca < v < ce, v, ?)]], SHADE(u , v, x, a, b, n)) 


STR(xl , yl , x2 , y2 , 1 , p) creates a broken line connecting points Pi(xi ,ji) and P 2 (x 2 




"P" 


0RDNER(x , y , 1 , p) creates a pair of broken lines parallel to the axes containing P(x,jp): 


(*,y) 


X 
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In my last figure I used additionally LABEL.MTH from DNL#16: 



Figure 5 

f4(x) := 0.3- (x + 2.5) • (x + 1) • (x - 2.7) + 1 

2 

f5(x) := 0.05-C3 - 2-x)-(x + 4-x + 5) 

SHADE(f5(x) , f4(x) , x, -3.5, 2.8, 50) 
S0LVE(f4(x) = f5(x), x) 


x = -0.763 v x = -2.65 v x = 2.19 


[STR(-2.65, 0, -2.65, f4(-2.65), 0.1, 0.05), STR(-0.673, 0, -0.673, f4(-0.673), 
0.1, 0.05), STRC2.19, 0, 2.19, f4(2.19), 0.1, 0.05)] 


It is not very difficult to extend the given functions: so you could add scale factors to consider differ- 
ent scalings for the two axes in ORDNERQ . 



I found this nice snail-shell in a Mathcad book and I tried to produce it with DERIVE. 
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Visiting DERIVE 
in Austrian Schools 


Some Austrian teachers 

1 DERIVE in a 3 rd form (age 13), Alfred Eisler, Tulin 

Some impressions and results of using DERIVE in a 3 rd form of a gymnasium. 

Use of the program 

We don't use DERIVE continuously, but we use the PC if its usage is “fitting well 64 to the teaching situation. 
Knowledge of the pupils 

Round 30% of the pupils are using a PC at home, but none of them has DERIVE available. Until now we have 
used DERIVE in connection with two items: 


1. Solving equations 

2. Working with co-ordinates. 

to 1) Solving equations was used to introduce the program and to discuss the most important commands (Au- 
thor, Simplify, Solve and Remove) and to show some of DERIVE's capabilities. We talked about the dif- 
ference between working with DERIVE and working with a pocket calculator. 

to 2) The pupils had to find the co-ordinates of points, Author them and plot the resulting figure. Between 1) 
and 2) there was an interval of more than a week. 

I made the following experiences: 


• The pupils like to go into the PC-lab. 

• The soLve-command creates a very strong “Aha - effect 44 . 

• The biggest problems occured in editing the expressions. Pupils make many mistakes - they don't find the 
characters on the keyboard and that causes different expressions and then different results. 

• We need too much time for editing and then correcting the input. 

• The pupils are forgetting the most things within one week. I had to start explaining the first steps from the 
beginning in the next lesson. If we don't use the PC continuously then it would be better to work in 2 hours' 
lessons. 


The last three points seem to be the problem. Until now I cannot give an answer to the question “What is the PC 
in this age good for (in maths, of course)? 44 , which is satisfying me. I shall report again in the future. 


II & 

#1: "Beispiel 708 j Lehrbuch Laub" 

# 2 : [- 1 , -21 
#3: 11, -21 

#4: 11, 1] 

#5: 12j 1] 

#6: 10, 31 

#7: [-2 j 1] 

# 8 : [- 1 , 1 ] 

#3: "The points haue to be'' 

#10: "connected by hand." 

#11: ^ 
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2 With DERIVE on the floor (age 12), Bernhard Wadsack, Vienna 

DERIVE and large prime numbers 

In connection with the chapter prime numbers in teaching mathematics in a 2nd form the pupils mostly 
are very astonished that there is an infinite number of prime numbers. The fact that there exists no 
largest prime number and that you are able to proof this thesis is nearly incredible for them. 

At the other hand for every moment in maths' history exists a largest known - and published - prime 
number. A special part in hunting larger and larger prime numbers were playing the so called 
Mersenne prime numbers of the form 2 P - 1 withp e {2,3,5,13,17,19,31,61,89,107,127,521,607,1279, 
2203,2281, ....}. Using computers additional prime numbers could be found among the Mersenne 
numbers. In 1986 it could be shown, that 2 P - 1 with p = 216091 is a prime number. This number has 
65050 digits. 

Now I tried whether DERIVE 2.56XM is able to calculate this huge number. First I quickly reduced 
the decimal places to 2 to harass DERIVE not too much and then I was curious whether DERIVE 

would Simplify 2 216091 - 1.1 waited and waited, and after 43.5 sec it was done. DERIVE had 

mastered the calculation. Now it took some time until an endless line appeared on the screen. There it 
was at last, the - until 1986 - largest known prime number. 31 pages - thirty one in letters - full of dig- 
its, that were unfolded 9m of printer output. In that very moment I had the idea to expand this huge 
prime number in the classroom of my 2nd form. 

The pupils made big eyes when I told them about this enormous prime number, but their astonishment 
became even greater, when I started to unfold the printout on the floor. To provoke a longer activity 
with this number I set them the task to look for chains of same digits or sequences of digits like 123.... 
All the pupils dispersed at both sides of this “number carpet 44 and doggedly searched for columns of 
same digits: soon we heard messages of success about groups of three digits the same, then sporadi- 
cally about groups containing 4 digits and it lasted a while until one pupil had discovered five fours. 

You easily can imagine that this lesson ended too soon for many pupils. Many of them had liked to 
search a longer time when I had to pick up 44 my number 44 in the break. The lesson had been very inter- 
esting, the pupils were enthusiastic. Whenever had it been, that they could glide on the floor during a 
maths lesson. 

I am sure that they will not forget this lesson for a long while. In the next lessons many questions fol- 
lowed, e.g. whether they themselves would be allowed to work with this program which is able to 
handle with such large numbers. I will comply their wish with pleasure because I've intended to use 
DERIVE in maths teaching since long. 

Comments: You can Simplify 2 216091 - 1 only with DERIVE XM. (I worked on a 386, 8 MB RAM). I 
couldn't calculate 2 500000 in Exact Mode. 

-Sincerely Bernhard W. 

(Enclosed were 31 pages ~9m ~ 120g: 74 6093 2844 1, Josef) 

(With DERIVE 6.10 2 216091 -1 takes 0.64 sec and 2 500000 takes 3.41 sec = 40 pages 

On the next two pages you can find two working sheets which I use to get the pupils accus- 
tomed with DERIVE and to develop some maths skills. On the first paper they hopefully 
achieve a feeling for the structure of the expression and the different levels of the subex- 
pressions. In the second paper they can see the expressions on the screen - they see and 
they learn to see the mistakes they have made in setting wrong parentheses. When they 
later will work with their pocket calculator they will not see the expressions, they often fail in 
calculating similar numerical expressions. Additionally they have to make a difference be- 
tween decimal places and significant places. 
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3 Rediscovering structures with DERIVE (age 16 ), Josef Bdhm, Wiirmla 

DERIVE/lll/1 .2 

a 2 +2b 2 3a+4b ^2a+3b j 

(a+3b) 3 5b + c 2 3a 3 +(a-4b) 4 

Edit this terrible rational expression! Compare carefully the result with the expression on the paper. 

Then highlight using the arrow keys the designated subexpressions. Switch off the Line Edit Mode in 
the message line (using the F6 - key). Start always with . i will bring you one logic level deeper 
into the expression; with — » and <— you can move within this level. Using T you will climb up one 
level. (You can help yourself with T, T, at any position.) 

Note how many keystrokes you needed: 

a 2 +2b 2 3a+4b (2a+3b j 

(a+3b) 3 5b +c 2 3a 3 +(a-4b) 4 

a 2 +2b 2 3a+4b ^2a+3b j 

(a+3b) 3 5b+c 2 3a 3 +(a-4b) 4 

a 2 +2b 2 I 3a+4b (2a + 3b 2 ) 

(a+3b) 3 5b+c 2 3a 3 +(a-4b) 4 

a 2 +2b 2 3a+4b ^2a+3b j 

(a+3b) 3 5b+c 2 3a 3 +(a-4b) 4 

a 2 +2b 2 3a+4b (2a+3b j 

(ia+3b) 3 5b+c 2 3a 3 +(a-4b) 4 

3+2b 2 3a+4b + (2a + 3b 2 ) 

(a+3b) 3 5b+c 2 3a 3 +(a-4b) 4 


a 2 +2b 2 

3a+4b 


(2a+3b 2 ) 3 

(a+3b) 3 

5b +c 2 

1 

3a 3 +(a-4b) 4 

a 2 +2b 2 

3a+4b 


(2a + 3lP) 3 

(a+3b) 3 

5b +c 

1 

3a 3 +(a-4b) 4 

a 2 +2b 2 

3a+4b 


(2a + 3b 2 ) 3 

(a+3b) 3 

5b +c 2 

\ 

3a 3 +(a— [^b) 4 
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DERI VE/II 1/1.3 

Calculate the numerical expressions using your pocket calculator. Note the precision wanted: 
(n d = n decimal places; n s = n significant places). Check your results using DERIVE. You 
have to use the approX command. Pay attention to Option Precision and Notation and the 
Digits. 



Note: The n-th root expression must be entered in DERIVE as (expression) A (l/n). 
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4 Dynamic Systems (age 18), J. Lechner & H. Achleitner, Amstetten 

Josef Lechner and Helmut Achleitner are teachers at the gymnasium in Amstetten, Lower Austria. 
This is one of Austria's “DERIVE show - schools". They both are working since a long time using 
DERIVE. Because this school is an “Informatics Gymnasium “ they spend more time for mathematics 
tuition and they also can give more attention to the programming aspect. Supported by their teachers 
the students develop among others tools with DERIVE, which they later on can use as a “Black Box“. 
See on this page one of these tools concerning dynamic systems. 

Josef & Josef have sent several examples of examination papers. I tried to translate them. So, if there 
are any mistakes, then they are my fault, don't blame Josef & Helmut. (But blame Josef _2, that’s me.) 

I didn't translate the toolbox for the Dynamic Systems, „Wachstumsprozesse Vernetzte Systeme “ 
(Growth processes, Cross-linked (= dynamic) systems) 


Wachstumsprozesse 
Vernetzte Systeme 

Wertetabelle fur explizite Funktionen 

• F(x) := Funktionsterm 

• TAB(xO , sw, n) := VECTOR([k, F(k)], k, xO, xO + (n - l)sw, sw) 

Wertetabelle fur rekursive Funktionen 

• F(x) : = Rekursionsgleichung 

• TAB_R(xO , n) := VECTOR([k, ITERATE (F(v) , v, xO, k)], k, n) 

Vernetzte Systeme (Phasendiagramm) 

• F(sys_parl, sys_par2) := [System 1, System 2] 

• TAB_PH(st_l , st_2 , n) := ITERATES(F(ELEMENT(v, 1), ELEMENT(v, 2)), 

v, [ st_l , st_2] , n) 

Wertetabelle fur die Systemvariable 1 

• TAB_PARSYSl(st_l , st_2, n) := [VECT0R(k, k, 0, n) , 

ELEMENT(TAB_PH(st_l , st_2 , n) ' , 1)]' 

Wertetabelle fur die Systemvariable 2 

• TAB_PARSYS2(st_l, st_2, n) := [VECTOR(k, k, 0, n) , 

ELEMENT(TAB_PH(st_l , st_2, n) ’ , 2)]' 

Darstellung im x„ - x n+ i - Diagramm 

• TABS(xO , n) := ITERATES ([ELEMENT(v, 2), F(ELEMENT(v, 1))], v, 

[xO, xO], n) 
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The examples: 

1 . Cats and Mice 

A region is populated by cats and mice. 

a) The “Standard Moder 

The number of mice is determined by its natural increase of 30% annually and its decrease caused 
by their enemies, the cats. This decrease is proportional to both the number of cats and mice with a 
factor of proportionality b = 0.005. 

Because each of the cats claims a certain district only 150 cats can survive in the whole region. 
Hence the increase of the cats is proportional to the number of mice and to the existing growth po- 
tential (proportionality factor d = 0.0001). The cats' natural annual decrease is 15%. When obser- 
vation starts there are living 2000 mice and 40 cats in this region. 

State the difference equations and investigate the long-time evolution of the mice' and cats' popu- 
lation. Will there appear any position of balance? If yes, then calculate it and explain the reason for 
the existence of equilibrium. 

(Plot the graph and save as CAT1 !) 

b) “Immigration of Cats 44 

In the region observed there are a lot of cat lovers, so annually z animals are added. At the other 
hand it can be assumed, that the cats are loosing partially their natural instinct with the result, that 
the increase of the “Standard Model 44 is only proportional to the number of the mice and cats (and 
does not longer depend on the growth potential). 

Form again the difference equations and investigate the behaviour of the system (take z = 3 and 
don't change the other parameters). What can be observed? What is the effect of an increasing z? 

(Save as CAT2.) 


2. A recursive sequence 

Show: The sequence <x n > defined by the recursive equation (difference equation): 


1 


X +1 2 


x + 


a > 0 


X 

n / 


with x 0 as starting value converges to x = Va . 

Represent the sequence in a x n - x n +i - diagram. (x 0 = 10, a = 2). 
(Plot and save as HERON1.) 


Explain why this difference equation is not linear. 






D 
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3. Life in the Soil 

Life in the soil causes an annual increase of the concentration of nutritive substances of 32% (i.e. 
forming humus). Grain cultivation consumes these substances so without fertilizing the ground less 
grain is growing and the extraction of the nutrients - starting with a 400 units quantity - will show a 
decrease rate of 5% annually. 

a) How will the quantity of nutrients x n develop within the next years, if there are now 1070 units of 
quantity in the ground? Could it be possible to exploit totally the soil, i.e. the nutrients would de- 
crease to 0? (if so, then find the time.) 

b) How many units of quantity could be consumed annually by cultivating grain that the quality of the 
soil could remain the same? Give reasons! 

c) What is the type of the difference equation presented and calculate the first four elements of the 
sequence. 

d) Annual increase of nutrients is not constant 32% but will reduce annually by 1% (31%, 30%, ....). 
Set up the corresponding difference equations. 

(Save as GROUND 1) 

The three examples from above are accomplished with a fourth one from probability theory and form a 
2 h - written examination (100 min) in a 7th form (age 17). Here are two more examples, one of them 
dealing with financial mathematics: 

4. Fish and Fishermen 

A population of 950 fish has an annual growth rate of 50%. The annual catch numbers are increasing 
by 10%, starting with 400 fish in the first year (then 440, 484, ). 

a) Describe the growing process using a difference equation. 

(pc n is the number of fish at the end of year n ; x 0 = 950). 

b) How many fish will represent the population in the 3rd, 4th and 5th year? 

c) Using DERIVE simulate the process for the first ten years. 

d) When will the population die out? Reduce the catch rate, so that the fish will not die out. 


5. Mr Akiar and his loan 

Mr Akiar takes a loan of AS 250000 with an interest rate of 1 1% annually. Exactly one year later he 
starts repayment. His first payment is AS 25000. His trade union has promised an annual 8% increase 
of salary, so Mr Akiar has the intention to raise accordingly his repayments also by 8% annually. 

a) Give a representation of the debts in respect to time (as a value table). Will the bank director accept 
Mr Akiar's offer? Will Mr Akiar be able to finish repayment? Write down reasons for your an- 
swers! 

b) Use DERIVE to represent the development of Mr Akiar's debt. 

c) Union's negotiations are failing year for year and the promised rises of salary don't become true. 

At how many percents has the annuity at least to be increased annually that the loan can be repaid? 
What is the “limit interest rate“, that repayment is yet possible (2 decimal places). 
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At last I'd like to present an example from an end examination incl. solution. In Austria we call this end 
examination “Matura“, what is the “Abitur“ in Germany. When a student has passed the Matura, then 
he/she is “maturus“ (I at., = ripe) to attend university. 

6. Growth under self poisoning (Environment pollution) 

A biotope has food for round 20000 individuals maximum. This population produces metabolic prod- 
ucts, which are leading to a decrease of the population. The natural growth rate is proportional to the 
existing growth potential, the rate of decrease caused by the metabolites is proportional to the number 
of the present individuals and to the amount of the metabolites. 

The increase rate of the metabolic products is proportional to the number of individuals. These prod- 
ucts are decomposed by bacteria. The decrease rate is proportional to the existing amount of the prod- 
ucts. 


x n number of individuals at time n (x 0 = 5000) 

s n amount of metabolites at time n (s 0 = 0) 

z increase rate of the individuals (proportionality factor of the 

growth potential) (z = 0.3) 

a decrease rate caused by the influence of the toxic metabolites (a = 0.001) 
c increase rate of the metabolites (c = 0.01) 

d decrease rate of the metabolites caused by the bacteria = 0.1) 

Describe the system using difference equations. Investigate the long-time behaviour of the popula- 
tion's quantity and the amount of the metabolites. (Give the first 50 pairs of values!) 

Give an interpretation of the result. 

Is there any equilibrium possible? Make your conclusions from the long-time-simulation. 

Express the fix points by the parameters. 

Give a graphical representation using an appropriate tool. 

Give a simulation for the special case, when d = 0 and interpret the outcome. 

Comment: Save the simulation on a diskette. Produce a hardcopy of the simulation and the graph. 
Label the plot by hand. 

Solution: 

The model : x n+ \ = x n + z (20000 -x n )- a • x n - s n 

Sf 2 +i Sn c x n d • s n 


#1 : L0AD(D : DYN_SYS . MTH) 

#2: F(x, y) := [x + z-(20000 - x) - a-x-y, y + c-x - d-y] 

#3: [z := 0.3, a := 0.001, c := 0.01, d := 0.1] 


#4: 


TAB_PH(5000 , 0, 50) 
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5000 

9500 


0 

50 


#5: 


4 

1.2175-10 

6390.046754 

6390.013534 

6389.986451 


140 

638.9654646 

638.9693857 

638.9725824 



#6: TAB_PARSYS1(5000 , 0, 50) 

‘ 0 5000 

1 9500 


# 6 : 


4 

1.2175-10 


48 6390.046754 

49 6390.013534 

. 50 6389.986451 . 

#8: TAB_PARSYS2(5000 , 0, 50) 

(TAB_PARSYS2(5000 , 0, 50)) 

#9: [48 51] 


#10: 


47 638.9606551 ' 

48 638.9654646 

49 638.9693857 


50 638.9725824 J 
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Interpretation: After a big increase of the population at the beginning (black), we can see the conse- 
quences of self poisoning and after a while appears a state of equilibrium. 


The values of equilibrium : 


x n = x„ + z (20000 - x„) - a ■ x n ■ s n 
s„ = s n + c ■ x„- d ■ s„ 


Solve this system of equations for x„ = x and s n = s 

S0LVE(x = (F(x , y)) a y = (F(x, y)) , [x, y]) 

# 11 : 1 2 

#12: (x = -9389.87 a y = -938.99) v (x = 6389.87 a y = 638.99) 

Equilibrium in x = 6340 and y = 639 
Case d = 0 
#13: d := 0 

#14: TAB_PARSYS1(5000 , 0, 50) 

#15: TAB_PARSYS2 (5000 , 0, 50) 



It seems to converge to an equilibrium! 

S0LVE(x = (F(x , y)) a y = (F(x, y)) , [x, y]) 
#16 : 1 2 

#17: false 


But there is no solution for the equilibrium. Let's investigate the long time behaviour! 
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#18: TAB_PARSYS1(5000 , 0, 100) 

#19: TAB_PARSYS2(5000 , 0, 100) 



#22: TAB_PH(5000 , 0, 50) 



#22: TAB_PH(5000 , 0, 50) 
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IL 


J 


Treating the problem with the TI-92 / Voyage 200 




jfLDTS 


^ ul=ulCn - 1) + z-(2GGGG - ul(n - 1)) - a ulO 
ui 1=5080 


u2=u2(n ■ 
ui2=0 
u3= 
u i 3= 
u4= 
u i 4= 
u5= 

ui 5= 


l) + c- ul(n - 1} - d ■ u2(n - 1) 


ul < n > = ... -1 > > -attnl < n -1 > *li 2 < n -1 >1 

Mil IN FiliD EKACT SEG 







IT1M 



0.0000 




mi'li'MiiEt 

^uluWulnlulum 

i in" 































mm 

mmm 

|w'V4- 

*v~T 1 

In 

ul 

u2 



EMawwasM 



EEHHMSM 




145.0000 

6300. 10SS 

63S.0475 


146.0000 

6300. 1375 

63S.054S 



laaaBiiBEaBM 








ft3R0. ORA!=i 

WIWIM* 


ul<n>=6389 -9864512702 


MAIN 

RAD EKACT 

SEG 


We define the two sequences in the Y=Editor (first of all turn to the Sequence Mode!). The table gives 
the numbers of the individuals for the following years. Compare with the DERIVE results. 





jpLDTS 
^ ui=|^ 
uil = 
u2=U 
Ui2=q 
u3= 
u i 3= 
u4= 
u i 4= 
u5= 

ul<n>=ul<n-l> + B^<20000-ul<n-l.. 

USE * AND * TD DF'EH CHDICES 


hxes 

fixes CUSTOMS 

i-J 

X Axis. . . ul-> 

V Axis. . . u2-> 

£Ent-er-=5flUEl > CE5C=CAHCEL^ 


We can plot the two sequences presenting the number of individuals and metabolites. 



The plot of the phase diagram can be achieved by customizing the axes - either in points or as a line. 




This are the respective plots for the special case d = 0. 
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5 Let's work in pairs (age 18 ), Hubert. Voigt & Karl Ohlinger, Perg 

Hubert and Karl are both teachers at a “Handelsakademie“, the same type of school on which I am 
working as a teacher. So I am very interested how they are dealing with the phenomenon “Teaching 
maths with a PC“. They have made the experiment to allow their pupils to work in pairs in their as- 
sessments. And they think, that it was successful. The marks didn 't change significantly, but the stu- 
dents reaction was enthusiastic. Hubert and Karl could observe a considerable learning process in 
developing teamwork and cooperation. 

See one of their tests: 


3. Test. IVA, Group A, 20.4.1994 

Pay attention: you have to write down the results on your paper (try to do it in words, like “The profit 
is “). Short comments between the calculations are welcome. 


What you have to work through together: 

1 .a) f(x) = ax 3 + bx 2 + cx 

The graph shows a maximum for x = 0.5486, a zero at x = 1 and a point of inflexion at x = -1/3. 
Find the equation of the function. (y = x 3 + x 2 - 2x). 

b) The graph from a) is intersected by the graph of y= ^ jx + 6 - 2 . 

Calculate the area between the graphs. 

2. A tank must hold 500 litres. Its length is 4m and its cross section is a rectangle with a semicircle 
at the bottom. A cover plate is not needed. Which dimensions minimize the consumption of mate- 
rial? 

3. Given are a supply function p$ = 0.25 x 2 + 3 and a demand function p D = 10 - e 0 4x . 

a) Give a sketch of both functions, according to the plot on the PC-screen. 

b) What is the new position of the equilibrium if there is a 10% subsidy? 

Are there any surplus goods? What is the reason for that fact? 

c) Which tax rate is necessary to give a total tax yield of 5 monetary units? 


This part has to be worked through separated (you may help each other). 

4. a) Find the price elasticity for the demand function p D (x) = 10 - ln(4x). (without using the PC) 
Simplify the result. 

b) What is the elasticity s D for x = 5 quantity units. Interpret the result for a 10% increase in price. 
What are the new price and the new demand? Check s D using the new data. 

c) Explain the conception ’’elasticity”. When can we speak of an ’’elastic” demand? 


5. 


Find the antiderivative of 



Write down the result using root signs. 


The partners were: 


I believe that example 3 could be interesting for those of you who are not accustomed with 
economic mathematical models. So I asked Hubert to send his solution: 
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Solution for 3.): 

The intersection point of the two curves is (x « 3.46,/? « 6.00). That is the equilibrium point 


10% Subsidies: 


Subsidies might lower the supply-price, if they are transmitted to the consumers, so with 

p s (x ) = (0.2 5 X 2 + 3) • 0.9 we obtain the new equilibrium: 

(0.25x 2 + 3) • 0.9 = 10 - e°' 4x => x = 3.65 QU and/? = 5.69 MU 

If subsidies would raise the producent's income, then they lead to an increase of the quantity of sup- 
ply. 

The new price is 6.00 (from the equilibrium) • 1.10 = 6.60. We find the new quantity of supply: 

6.60 = 0.25x 2 + 3 => x = 3.80 QU. 

The sold quantity is 3.46 QU, so we obtain a difference of 0.34 QU. (Surplus of grain, milk, etc.) 

Total tax yield of 5MU: 

Taxes raise the supply price. /? is the tax rate 

Ps (x) (1 + /?) =/?d(x) => /? s (x) • p =p D (x ) -/? s (x) = the tax/QU (*) 

Total tax = quantity • tax/QU => 5 = x • (/? D (x) - /?s(x)). Solutions: xi = 0.94; (x 2 = 2.89). 

Substituting for x in (*) and solving for p yields: p\ = 1.65, i.e. a taxation rate of 165% and the 
second solution p 2 = 0.34 = 34% 


Together with some other tests Hubert sent a paper as a result of an idea from DNL#13 (The 
circumscribed circle of a given triangle in CAS and Spreadsheets, D.Sjostrand). Hubert pre- 
pared two teaching units - one for the circumscribed, the other for the inscribed circle. Here 
is the second one: 


The inscribed circle of a triangle 

The calculation of the incircle is more difficult. You will need some information. 

The centre of the incircle is the intersection point of the angle bisectors. You will obtain its radius if 
you draw the perpendicular line from the centre with respect to one side. 

Task 1: Draw the incircle into the given triangle: 
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To find the equation of the angle bisectors you must know, that each angle bisector divides the op- 
posite side proportional to the adjacent sides. 


Task 2: 


In the triangle above measure the sides a, b, c; and the sections a ? , a”, b ? , b", c ? , c" and calculate the 
ratios. 

Example: c : b = a ? : a". 



We obtain 6 equations for the co-ordinates 
of the 3 dividing points Ti >2 , 3 . But we need 
only 4 of them: 

(x 4 - xi) : (x 2 - x 4 ) = b : a 
(x 5 - x 2 ) : (x 3 - x 5 ) = c : b 
(x 6 - x 3 ) : (xi - x 6 ) = a : c 
(y 4 ~y\) : (y2-y4) = b : a 
Os -yi) : Cf3 -ys) = C : b 
(ye-y3) : (yi-y6) = a : C 


The solutions of this system of equations are the co-ordinates of the dividing points Ti,T 2 ,T 3 . 
Note that these points are not the osculating points! ! 


Task 3: 

Given is the triangle A(-5/-3), B(5/-4),C(2/4). Draw this triangle on a separate sheet of paper. Calcu- 
late the point T 3 (x 4 /y 4 ) and add it to your sketch. Notice: the lengths of a, b and c are the distances of 
two points. 

Calculation of the other points is not a big problem, but now we will use DERIVE. For this purpose 
load a so called Utility file by File > Load > Utility File: i nci rcl e . mth. You cannot see the file and 
the formulae. 

You have to Author the co-ordinates of the vertices: [xl := - 5, y 1 := - 3, x2 := 5, etc ]. 

The following functions are provided for you. 

bs_a bisector of a 

bs_b bisector of p 

bs_c bisector of y 

radi us the radius of the incircle 

centre the centre of the incircle 

ci rcle the incircle 

Author the name of the function needed, then Simplify or Approximate. 

Task 4: 

We want to calculate the angle bisectors - remember the equation of a line which passes through two 
given points (2nd form). 



Find the equation of one angle bisector using a vertex and the opposite dividing point: 
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Then let DERIVE calculate all the angle bisectors and check your result(s). (Use Expand) 

Task 5: 

Take any two of the bisectors and calculate the intersection point. Compare your result with 
DERIVE 's answer on the provided function CENTRE! 

At last plot the triangle, the dividing points, the bisectors, the centre and the incircle. For plotting the 
triangle set Option > Display > Points > Connected > Yes. Then set this option back to Discrete. 

If you like you can see the "DERIVE - program”, with Transfer Load Derive: incircle. 

#1: inci rcle.mth 

#2: InputMode := Word 

#3: [xl := -3, yl := -2, x2 := 5, y2 := -1, x3 := 2, y3 := 6] 

' xl yl 
x2 y2 

#4: triangle := 

x3 y3 
. xl yl 


2 2 


#5: 

a := 

7((x 2 - x3) 

+ Cy2 - y3) ) 





2 

2 



#6: 

b := 

VCCxl - x3) 

+ (yl - y3) ) 





2 

2 



#7: 

c := 

V((xl - x2) 

+ (yl - y2) ) 



#8: 






#9: 

Radi us : 






a + b + c 




#10: 

s := 







2 






V(s-(s 

- a) • (s - b) • (s - c)) 



#11: 

radius := 







s 






r f x 4 - xl 

b 

\\ "N 


x4 : 

- IF xl = x2, 

xl, SOLUTIONS 

— J 

x4 

#12: 



L L x2 - x4 

a 

JJ 



V 



1 ) 



r 

r r y4 - yl 

b 



y4 := IF yl = y2, yl, SOLUTIONS = , y4 

#13: l l y2 - y4 a 



x5 := IF x2 = x3 


SOLUTIONS 


x5 - x2 





#14: 


x3 - x5 


1 ) 






D 
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y5 - y2 c Y 

\ A 

y5 := IF 

y2 = y3, y2 , 

SOLUTIONS 

= — , y5 


#15: 



y3 - y5 b ). 

/ 





lj 




s 

x6 - x3 a Y 

\ \ 

x6 := IF 

xl = x3, xl, 

SOLUTIONS 

X 

II 


#16: 



s. xl - x6 c >. 

/ 





i) 




s 

y6 - y3 a Y 

\ N 

LL 
1— 1 

1 1 

<X> 

> 

yl = y3, yl, 

SOLUTIONS 

VD 

> 

II 


#17: 



yl - y6 c ). 

/ 





lj 



‘ x4 

y4 ' 

#18: d_points := 

x5 

ys 


. x6 

y6 . 


#19: 


# 20 : 


# 21 : 


# 22 : 


bs_a := IF 


bs_b := IF 


bs_c := IF 


y5 - yl 

xl = x5, x = xl, • (x 


x2 = x6, x = x2 , 


x4 = xB, x = xB, 


x5 - xl 
y6 - y2 
x6 - x2 
y4 - yB 


(X 


(x 


x4 - x3 

m := (S0LUTI0NS(bs_a = bs_b, x)) 


- xl) + yl 

- x2) + y2 

- x3) + y3 


n := lim bs_a 
#23: x->m 

#24: centre := [m, n] 


2 2 2 
#25: circle := (x - m) + (y - n) = radius 


Note: Vertical bisectors can be plotted easily in 
DERIVE 3.x because of its implicit plotting capa- 
bility. In earlier versions you can help yourself by 
creating a special plot_bisector function: 
pbs_a:=[[x1 ,y1],[x5,y5]], etc. In this file you find 
the circle given in implicit form, too. 
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Titbits from Algebra and Number Theory 

by Johann Wiesenbauer, Vienna 

As this is the first article of my series after the arrival of DERIVE 3, 1 am going to try out some of the new fea- 
tures of it. My first impression is that the performance of DERIVE has been enhanced greatly by a lot of very 
useful functions, even though the first official versions of DERIVE 3 contained some tantalising bugs (as a 
matter of fact, more serious ones than the Beta version 3.0y). In particular, some of the following routines will 
only work properly, if you have version 3.04 of DERIVE or later. 

Speaking of bugs, the DERIVE-routine for the ju-function in DNL #15 yields the wrong value 0 for n=l and 
hence should be changed to 

p(n) := IF(n = 1, 1, IF(DIM(FACTORS(n)) = X( ( FACTORS (n)) u2) , 

(-1) A DIM( FACTORS (n)) , 0)) 

The reason for this is a tiny inconsistence in the current DERIVE: factors(factor(l)) returns [1] instead of the 
proper value [ ] . 

The following simple function derived from p is also sometimes useful: 

SQUAREFREE(n) := IFCn(n) A 2 = 1, true, 0) 

Let us turn to more demanding tasks now, e.g. what about the DERIVE-routine for the number 8(n) of divisors 
of n?" Try it on your own first! 

6(n) := DIM(SORT(VECTOR(nCu_) , u_, { [1] } •n(VECT0R(MAP_LIST([v_4,l A k_] , k_, 

{0 v_j,2}) , v_, FACTORS CrO))))) 

VECT0R(5(i ) , i , 20) = [1, 2, 2, 3, 2, 4, 2, 4, 3, 4, 2, 6, 2, 4, 4, 5, 2, 6, 2, 6] 

One of the most important functions of number theory is certainly Euler's cp-function which can be programmed 
as follows: 

()>(n) := n-[T(l - l/v_j.l, v_, FACTORS(n)) 

VECTOR ((|>(i ) i i, 20) 

[1, 1, 2, 2, 4, 2, 6, 4, 6, 4, 10, 4, 12, 6, 8, 8, 16, 6, 18, 8] 

(Note that because of the same DERIVE-bug as above the case n=l needs an extra definition.) 

Among others, the cp-function can be used to calculate the inverse of a mod m, if m is of moderate size: 

INV(a, m) := IF(GCD(a, m) = 1, M0D(a A (c|>(m) - D , m)) 

INV(2 , 1234567891234567891234567) = 617283945617283945617284 
M0D(2- 617283945617283945617284, 1234567891234567891234567) = 1 

Of course, the inv-function based on Euclid's algorithm is far superior to the preceding one (0.1 sec versus 6.3 
sec for the example above!) : 

INV_(a, m) := IF(GCD(a, m) = 1, M0D( (ITERATE (IF(M0D(a_, b_) = 0, 

[a_, b_, c_, d_] , [b_, M0D(a_, b_) , d_, c_ - FL00R(a_, b_)-d_]), 

[a_, b_, c_, d_] , [a, m, 1, 0]))f4, m)) 

Even so, in view of the huge power that had to be calculated, the first implementation worked surprisingly fast. 
The reason for this: My prayers as to the DERIVE-implementation of the powermod-function (cf. #14, p.29) 

have been heard at last! This must have happened very recently , judging by two facts: Firstly that ridiculous 

powermod-function in the DERIVE-library is still there, and secondly, what is more, they forgot to change the 
mods-function, too! But don't look a gift horse in the mouth, as the saying goes, I am sure that this will be reme- 
died before long. In the meantime take this version of mods: 
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MYM0DS(a, n, m) := M0DS(M0D(a A n , m) , m) 

Here is still another function, that makes use of the improved version of mod: 

LEGENDRECa, p) := IF(S PRIME(p), "Input error!", MODS(MOD(a A ((p - l)/2), p) , p)) 

Of course, Legend re (a, p) is the well known Legendre symbol (a/p) from the theory of quadratic residues, 
which is defined for all integers a and primes p. There is an important generalisation called after Jacobi which 
cannot so easily be implemented. Here is my solution: 

JAC0BI(a, b) := IF(MOD(b, 2) = 0 v b < 0, "Input error!", IF(GCD(a, b) > 1, 0, 
IF(a = 1, 1, IF(ABSCa) > b, JAC0BI(M0DS(a, b) , b) , 

IF(a < 0, (-l) A ((b - l)/2). JACOBI (-a, b) , 

IF(M0D(a, 2) = 0, (-l) A ((b A 2 - l)/8) • JAC0BI(a/2 , b) , 

IF(M0D(a, 4) = 1 v M0D(b , 4) = 1, JAC0BI(b, a), 

- JACOBI (b , a)))))))) 

As to the performance of these two functions again an example: 

LEGENDRE(10 a 60 + 1, 10 a 68 + 99) = 1 (0.015 sec) 

JACOBI (10 a 60 + 1, 10 a 68 + 99)= 1 (0.000 sec) 

Both results show that the congruence x 2 = 10 60 + 1 mod (10 68 + 99) is solvable. One of its solutions, which 
differ by the sign only, is given by 

c : =MOD((10 a 60+1) a ( (10 A 68+100) /4) , 10 A 68+99) 
as can be seen by 
M0D(c A 2 , 10 a 68+99)= 

1000000000000000000000000000000000000000000000000000000000001 

The general proof for this is so simple that I cannot resist the temptation to insert it here: If CL = 6 2 mod p for 


integers a,b and a prime p , which is of the form 4&+3 and does not divide a, then a A ((p+ 1)/4) (in DERIVE- 
notation) is always a solution of square root of a because of 



bP +x =b 2 -b p ~ x 


a mod p 


where we made only use of Fermat's Theorem. The case p = 4 k+ 1 is far more difficult and maybe I will talk 
about it another time. 

As a final example for an application of the new (power)modfunction in DERIVE, I am going to deal with 
Pepin's primality test for Fermat numbers F = 2 2m + 1 (m > 0) which states that F is prime for m> 0 iff (if 
and only if!) 

3 (F -- l)/2 = -1 mod F . 

And here the DERIVE-implementation of it: 

MYMODS(a,n,m) :=M0DS(M0D(a A n,m) ,m) 

PEPINTEST(m) : =IF(m=0 , true , IF(MYM0DS(3 , 2 A (2 A m-l) ,2 A (2 A m)+l) = -1, true, false)) 

VECTOR(PEPINTEST(i) ,i , 0 , 10) =[ true , true , true , true , true , fal se , fal se , fal se , fal se , 

fal se , fal se] 

The calculation above for the first 11 Fermat numbers took only 2.8 (!) seconds, even though the last of the 
tested Fermat numbers has already 309 digits! Let me conclude with this most impressive example for the per- 
formance of the new DERIVE. If you have any suggestions as regards possible improvements of the routines in 
this article or interesting topics I would be very pleased to hear from you! (jwiesenb@email.tuwien.ac.at) 
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G P Speck, Wanganui, New Zealand 

I am using DERIVE Version 2.6 an a IBM 386 compatible. Having just recently joined the DUG, I am 
not conversant with the material in the first 14 Newsletters. However, I offer the following two Notes 
which I feel may be of interest to the User Group if the comments in these notes have not already been 
covered in previous Newsletters. 

Note 1 : A fault in the DERIVE Precision Option? 

When using Options Precision Approximate Digits: 5 I find that on returning to a view of the Digits 
entry it has changed from 5 to 6. The same is true if I enter Digits: 11; it changes from 11 to 12; simi- 
larly, if 

Digits: 17 is entered the 17 changes to 1; in general, if 5 + 6n is entered for the Digits, then the 5n + 6 
changes to 6 + 6n. Either 5 + 6n or 6 + 6n significant figures may in fact be displayed for each of the 
entries 5 + 6n and 6 + 6n depending on the intermediate entries and the order of the entries 5 + 6n and 
6 + 6n. If this is indeed a fault, it is a VERY SERIOUS FAULT! Have any other readers encountered 
this difficulty? 

DNL: This fault does not occur in versions 3.x any longer. But thanks for your useful com- 
ment. I am sure that there are many of our members who are still using earlier DERIVE ver- 
sions. And many thanks for your contribution: Muller's Method. Your second note might be a 
nice addendum to a DNL#1 6 's article. 

Note 2: A note on plotting Periodic Extensions 

Students are usually intrigued by the following elementary result when it is first encountered. 

If function F is defined for all real x then the periodic extension on the whole real line of F restricted 
to the half-open interval [a,b) is given by 

F (x - (b-a) FLOOR ( (x - a) / (b - a) ) ) . 

A DERIVE example of this result is shown below. 


y 

2 

1 

y 

2 

-2 -i 

i 2 

-l 

-2 

-2 -1 

1 2 

-1 
■ -2 

#1: := x 2 , a := 0.25, b := 1^ 

#2 : FH ( x > := IF ( IM (x > = O AND 0.15 < x < 1, F<x>> 

#3 : FE(x) := F[x - (b - a > -FLOOR [— ^ — ]] 

#4 : mi 

#5: ' 'FR (x ) is plotted 1 ef t above - FE(x> is plotted right above" 
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H. Voiqt, Perq, Austria 

I. The cubic 30x 4 - 91x 3 + 15x 2 - 91x + 30 could be solved using DERIVE 2.60 in both modes Exact 
and Mixed as well. But now with DERIVE 3.x it does not work in Mixed Mode. The system will 
hang. 

2. Students tried to plot y(x) = x + ln(x-3) - ln(x-5), and they could see a plot in a region, where the 
function should not exist, because it is defined for real x > 5. Do you have an advice? 

x — 3 

DNL: Yes. I have. Obviously DERIVE simplifies the logarithms internally to in So it is 

x-5 

very informative for the students to define the domain - as they have to do, when sketching 

the graph by hand. 


U(x) : =x+LN (x-3 ) -LN (x-5 ) 
DEF (x) : =x>3 AND x>5 
F (x) : =IF (DEF (x) ,U(x) , ?) 
F (x) 



Dr.H.J.Kavser, Dusseldorf, Germany 

Eine Alternative zu IN VN ORM.MTH aus DNL#15 

Function for the inverse normal distribution 

INVNORM in DNL#15 


#1: INVN0RM(p, p, ct) := ITERATE 


a 


SIGN(ct)-ERF| 


72 - a 72-p "I 


2-cr 


2-ct ) 


2 2 2 2 2 
a /( 2 -cj ) + a-p/ci - p /(2 -ct ) 


, a, p, 5 


72 -e 




2-7 M 
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NORMAL-INV, the Challenger 

Definition of the function ^-inverse (OlNV) 

#2: [aO := 2.51551, al := 0.802853, a2 := 0.010328, bl := 1.43278, b2 := 0.189269, 

b3 := 0.001308] 

#3: W(y) := 2-LN(y)) 

(a2*W(y) + al)*W(y) + aO 

#4: V(y) := W(y) 

C(b3-W(y) + b2)-W(y) + bl)-W(y) + 1 

<t>INV(y) := 

Ify>0Ay<l 
#5: If y > 0 a y < 0.5 

V(y) 

- V(1 - y) 

error-estimating: |6 1 < 4.5-10 a -4 

Definition of function NORMAL INV 


#6: NORMAL_INV(y , n. cr) := |J + a-O)INV(y) 

Testl (accuracy) 

#7: VECTORCEy, INVNORM(NORMAL(y , 105, 7.5), 105, 7.5), NORMAL_INV(NORMAL(y , 


105, 7 

.5), 

105, 

> 

l 1 

A \ 
LO 

120, 127.5, 1.5) 

120 

119. 

.999 

120. 

O 

O 

UJ 


121.5 

121. 

.495 

121 . 

.503 


123 

122. 

.964 

123. 

.003 


124.5 

124. 

.329 

124. 

.503 


126 

125. 

.472 

126. 

.002 


127.5 

126. 

.309 

127. 

.502 . 


Test 2 (velocity) 


#9: VECTORCEy, INVNORM(NORMAL(y , 105, 7.5), 105, 7.5)], y, 120, 127.5, 1.5) 
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120 

119.999 ' 

121.5 

121.495 

123 

122.964 

124.5 

124.329 

126 

125.472 

127.5 

126.309 . 


0.031 seconds (2.3 seconds 1995) 

#11: VECT0R([y , N0RMAL_INV(N0RMAL(y , 105, 7.5), 105, 7.5)], y, 120, 127.5, 


1.5) 


120 

120.003 ' 

121.5 

121.503 

123 

123.003 

124.5 

124.503 

126 

126.002 

127.5 

127.502 . 


0.000 seconds (0.7 seconds 1995) 

#13: INVNORMCx, 0, 1) 

Takes about 30 seconds 

#14: NORMAL_INV(x , 0, 1) 

Plots in an instant. 




Plot of #13 


Plot of #14 
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Final Comment: 

Accuracy of function INVNORM could be increased by increasing the number of 
iteration steps but time of computation - in any case higher than applying NORMALJNV - 
would be increased simultaneously. 


Three Pictures from my 3D-Gallery 



Helix of a Rotoide 
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This is the DERIVE session: 

#1: Precision := Approximate 

#2: Notation := Decimal 

#3: NotationDigi ts := 3 

#4: The two parabolic cylinders 

#5: length: 5,5,2; grids: 50,50 

#6: I F C | x | < 1.5 v | y | < 1.5, IF(|y| > |x| a |x| < 1.5, - x-x, 

IF( | x | > | y | a | y | < 1.5, - yy, 0))) 

#7: A heli coni cal tube 

#8: r := 2 + COS(p) 

#9: HE(p , f) := [r.COS(f), r.SIN(f), SIN(p) + 0.8-f] 

7T 

#10: ISOMETRICS HE(p, f ) , f, , 6 -tt, 30, p, 0, 2 • tt , 30 

2 

' ( n 

#11: C0PR03ECTI0N ISOMETRICS HE(p, f ) , f, , 6 -tt, 30, p, 0, 2 • tt , 30 

^ 2 

#12: A helix of a rotoide - Rotoidenwendelflache 

#13: r_ := 3 - p-C0S(3.f) 

#14: RW(p, f) := [r_-C0S(f), r_-SIN(f), p-SIN(3-f)] 

#15: ISOMETRICSCRWCp, f ) , f, 0, 2 - tt , 60, p, -2, 2, 50) 

#16: COPRODECTION(ISOMETRICS(RW(p, f ) , f, 0, 2 - tt , 60, p, -2, 2, 50)) 

ISOMETRICS gives one family of parameter curves - plot with small connected points - and 
C0PR01ECTI0N gives the other family. 

The isometric projections of the surfaces in the 2D-Plot window had and still have their spe- 
cific charm - and investigating the various mappings have their mathematical value. 

At the other hand recent DERIVE version offers 3D plots of functions z = z(x,y) and of sur- 
faces given in parameter form - like the surfaces from above. 
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The two cylinders 



The tube and the helical rotoide (top view) 



This is the snail shell from page 30 as a 3D-plot 






